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ANALYSIS OF GEOMETRY AND DESIGN POINT 
PERFORMANCE OF AXIAL FLOW TURBINES 

PART I - DEVELOPMENT OF THE ANALYSIS METHOD 
AND THE LOSS COEFFICIENT CORRELATION 

by A. F. Carter, M. Platt, and F. K. Lenherr 

Northern Research and Engineering Corporation 


SUMMARY 


This report presents the development of a stream-filament analy- 
sis procedure and a correlation of tota 1 -pressure-loss coefficients which 
forms the basis for a computer program with which the geometry and design- 
point performance of axial turbines may be investigated. This report is 
the first part of a two-part report; the second part will present a com- 
plete description of the actual computer program. 

Since one of the principal features of the analysis procedure 
is the solution of the radial equilibrium equation taking into account 
radial gradients of enthalpy and entropy, the computer program based on 
this analysis will provide the turbine designer with the freedom to con- 
sider arbitrarily selected di s t r i but i ons of tangential velocities (or 
stator exit angles) and radial distributions of work output as analysis 
variables. In addition, wi th the incorporation of a total-pressure-loss 
coefficient correlation, the computer program which results from the analy- 
sis can be used for a systematic investigation of the performance of al- 
ternative turbine designs for specified design requirements. 



INTRODUCTION 


Until relatively recently, the majority of axial turbines has 
been designed using an extremely simple approach in which the flow at each 
of the design stations is assumed to have radially constant values of to- 
tal temperature and total pressure. With this assumption and the additional 
assumption that radial components of velocity are zero, the radial equilib- 
rium equation Is considerably simplified. One solution of the simple ra- 
dial equilibrium equation is the irrotational or free-vortex solution in 
which the axial component of velocity is constant, and the tangential com- 
ponent varies inversely with radius. The principal merits of the free- 
vortex design of turbines are that the design velocity triangles are read- 
ily obtained, and that acceptable levels of performance have been achieved 
when such designs are suitably bladed. 

The principal disadvantages of the free-vortex approach to tur- 
bine design are that the model is basically inaccurate and that the de- 
% 

sign method is very restrictive. Turbine test results from stages de- 
signed using the free-vortex method clearly show significant variations 
in axial velocity, total pressure, and total temperature. It is unrealis- 
tic to assume that there will not be significant radial variations of to- 
tal pressure and temperature when there are large variations in the aero- 
dynamic design parameters and the blade sections from hub to casing in 
typical turbine designs. The radial variation of the stage aerodynamics 
of free-vortex turbines is an illustration of the restrictive nature of 
the simple design approach. With a selected mean-line design of a given 
reaction and rotor deflection, the reaction and deflection of the remain- 
ing sections of the rotor blading are merely functions of the design radius 
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of these sections and, since no radial variation in temperature drop can 
be considered, the nondimens i ona 1 loading factor, a Tc/u 3- * , rapidly 

increases towards the hub section of the blade. Thus, the over-all geome- 
try of a turbine is frequently dictated by considerations of the hub load- 
i ng. 

Any free-vortex turbine can be considered as a series of simply 
related elements, and obviously the freedom of the turbine designer would 
be considerably increased if the dependence of the over-all design geome- 
try on a selected mean line could be reduced. Consideration of a number 
of elements over the radial extent of a design forms the basis for a 
stream-filament analysis of turbine designs; with this approach it is not 
necessary to select radially constant work extraction or restrict the 
variation of tangential velocities to a particular type. Unfortunately, 
a turbine design cannot be considered as a series of independent filaments 
defined by surfaces generated by the rotation of streamlines about the 
center line of the machine. The radial locations of the streamlines 
which define the individual filaments of the flow are dependent on radial 
equilibrium and mass flow continuity requirements at each of the design 
stations through the machine. Thus, an essential ingredient of the stream- 
filament analysis of a turbine design is the simultaneous solution of the 
radial equilibrium and continuity equations at each of the design stations 
through a turbine. However, once the decision is made to consider radial 
variations of enthalpy and entropy and to include the effects of meridio- 
nal components of streamline slope and curvature, no simple analytical so- 
lution of the radial equilibrium equation is possible. Hence, it becomes 
desirable to use a computer for the solution of the problem and it is, of 



course, one of the primary purposes of the present report to develop the 
analysis such that it forms the basis for a computer program. 

Even though the preparation of numerical procedures for the 
stream-filament analysis of a turbine design is of importance in itself, 
unless the procedure incorporates soundly based assumptions concerning the 
losses associated with the elements of the blading, the detailed calcula- 
tion of interblade row aerodynamics will be of little value. Hence, an 
essential part of the development of what is fundamentally a design analy- 
sis computing system is the development of a loss correlation which will 
be an integral part of the final computer program. Such a correlation has 
been developed as part of the over-all analysis. Since the basic require- 
ment is for the analysis of turbine design-point geometry, the correlation 
obtained is based on the interblade row flow angles and velocities and 
makes no reference to the detailed design of the blading. The detailed 
design of the blading will, of course, influence the over-all performance 
of a turbine, and the correlation of loss coefficient must be regarded as 
a datum level of loss which a turbine designer may factor as his experi- 
ence and knowledge of the probable blading dictates. 

While the model of streamline flow through a turbine represents 
a considerable advance from that used in the recent past, it is still a 
considerable simplification of the extremely complex flow which will occur 
in an actual turbine. For example, it assumes the flow is axisymmetric 
and ignores the secondary flow effects. Both are factors which make it 
unrealistic to assume that the flow which passes between selected adjacent 
streamlines at the turbine inlet will be identically that which passes be- 
tween these streamlines at some later design plane. While not attempting 
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to quantitatively assess the amount of mixing which occurs throughout a 
turbine, the analysis includes a method by which the effect of mixing can 
be simulated. An additional factor which Is considered in the analysis is 
the addition of coolant flows to the main stream. Here again, the model 
which is proposed is a relatively simple one which attempts only a first 
order correction to the over-all stream-filament approach to account for 
the addition of mass flow at a temperature level which may differ from 
that which exists in the main stream at the point of admission. No at- 
tempt is made to differentiate between cooling flow which results from 
disk cooling or blade cooling, be it transpiration or convective cooling. 

Report Arrangement 

The first section of the report outlines the stream-filament 
approach to turbine design and introduces the basic equations used in the 
analysis. Modifications to the simple stream-filament approach are then 
considered. These mod i f i cat i ons are undertaken to include in the analy- 
sis the effects of changes in specific heat through a turbine as the tem- 
perature level falls, the interfilament mixing which occurs within a tur- 
bine, and the addition of coolant flows in a high temperature application. 
The following section discusses the investigation which leads to the se- 
lected correlation of tota 1 -pressure- loss coefficient which is to be an 
integral part of the analysis program. In the f I na 1 sect i on , the analy- 
sis is developed to the point where It forms the basis of the computer 
program. 

The relationship between kinetic-energy-loss and tota 1 -pressure - 
loss coefficients and the coefficients of the differential equations to 
be solved are presented in the appendices. 



THE STREAM-FILAMENT APPROACH TO TURBINE DESIGN 


I ntroducti on 

The basic model used for the stream-filament analysis of a tur- 
bine design consists of a series of streamlines which trace the path of 
the flow from known conditions at the inlet to the first stage of the 
turbine to the final design station at exit from the last blade row. 
Making the assumption that the flow is axi symmetri c throughout the entire 
turbine, a series of streamlines can be selected at the machine inlet to 
define a series of annular streamtubes. With the position of the stream- 
lines defined so that adjacent streamline surfaces contain a known frac- 
tion of the total flow, the entire flow field throughout the turbine can 
then be considered as a number of annular elements. Since the flow is 
assumed axi symmetri c , the flow through the turbine can be represented by 
a meridional section as diagrammatical ly illustrated in Figure 1. 

In general, even for the simplified case in which axial symme- 

% 

try is assumed, the flow path through the stages of a turbine will be 
quite complex and will be influenced by a number of factors not the least 
of which will be the detailed design of the blading* However, for an an- 
alysis of the geometry and the design-point performance of possible tur- 
bine configurations, the principal objective is to define the blading re- 
quirement, and hence the analysis is restricted to design stations imme- 
diately upstream and downstream of the blade rows. 

With the turbine subdivided into a number of elements, the de- 
sign of the over-all turbine can be considered as that of a number of 
individual turbines. Hence, by applying the fundamental turbine design 
equations to individual sections of the turbine, the stream-filament 
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approach will permit the spec i f i cat i on of the design requirements to ele- 
ment of sections of the turbine analysis. These individual elements, how- 
ever, cannot be considered in isolation for these are interrelated by con- 
siderations of radial equilibrium and the over-all design requirements 
which include satisfying mass flow continuity at each design station and 
the specified total power output of each stage. Hence, the principal re- 
quirement of a stream-filament analysis of a turbine design is for the so- 
lution of the radial equilibrium equation; the solution of the design equa- 
tions in the streamline direction is a relatively simple task, which is 
merely an extension of standard turbine design practice to individual 
streamlines. Since the blade row elemental performance is likely to be 
a function of radius, and the work output may be specified to vary with 
radius, the radial equilibrium equation must take into account radial 
variations of entropy and enthalpy. In addition, as indicated diagram- 
matical ly in Figure 1, radial components of velocity and streamline cur- 
vature in the meridional plane will have to be included in the analysis 
of the flow field at each design station. 

In summary, for the s t ream-f i lament analysis, the turbine de- 
sign analysis proceeds along selected streamlines from known inlet condi- 
tions with the positions of the streamlines, other than those at the annu- 
lus walls, determined from continuity and radial equilibrium considerations. 
Although one of the objectives of the analysis will be to integrate these 
individual aspects, it is convenient to review each in turn. 

Analysis Along Streamlines 

The basic analysis performed along streamlines will differ little 
from that which is normally undertaken along the mean line of any axial 
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turbomachine. In fact, in only two respects can the design procedure be 
considered different from standard mean-line design practice. Firstly, 
the radial components of velocity are considered. Secondly, the specific 
heat at constant pressure, C/> , of the gas is assumed to vary from station 
to station through the turbine. The effect on the stat i on-to-s tat i on vari- 
ation of the specific heat is discussed in a later section of the report. 

Since the analysis is to form the basis for a computer program, 
it is convenient to divide the total flow into filaments having equal 
fractions of the total flow. An odd number of streamlines are selected 
so that the central streamline can be used as a mean streamline, which 
serves as the starting point for the solution of the flow field. This 
mean streamline serves essentially the same function as the mean line in 
conventional turbine design practice. The radial position of the stream- 
lines at each of the design stations and the streamline values of meridio- 
nal velocity will in practice be determined from radial equilibrium and 
continuity considerations. Initially only the boundary streamlines ( $ I 
and £T 0 f Figure 1) are known, for it is assumed that these follow the 
specified annulus contours. However, in the following discussion of the 
streamline equations used in the analysis, i t wi 1 1 be assumed that the 
radial position of all streamlines and their associated meridional ve- 
locities are all known quantities. 

For the purpose of the current analysis, which is performed in 
meridional and tangential planes, the absolute velocity vector is defined 
by two velocity components,^ and and two angles,^ and . The nomen- 
clature for the flow at a point Pina streamline surface in an axisym- 
metric flow in a turbine annulus is illustrated in Figure 2. In terms of 
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the axial, tangential, and radial velocity components, the analysis vari- 
ables are defined as follows: 


Meri di onal velocity, 


s# 

4- 

X 

> 

it 

(1) 

Absolute veloci ty, 

V = 

V v* + v£ 

(2) 

Flow angle, 

fi -- 

l/x 

(3) 

Streamline slope angle, 

A 

- Xa~Z! Va -/V k 

(4) 


Consequences of the introduction of radial components of flow 
into the analysis are that the conventional representat i on of the velocity 
triangles, which represent the flow on cylindrical surfaces, does not rep- 
resent the absolute velocities and that the values of station radius may 
vary from station to station. Sample velocity triangles for the central 
streamline (53) shown in Figure 1 are shown in Figure 3 to illustrate the 
nomenclature. I t wi 1 1 be seen that a superscript is used to denote rotor 
relative flow angles and velocities. Tangential components of velocity 
are assumed to be negative when opposed to the direction of rotation. 
Hence, rotor blade relative exit angles have negative values. 

As in any design, the calculations performed in the flow direc- 
tion are to determine the total pressure, total temperature, and the ve- 
locity components at each of the design planes; in the streamline analysis 
these quantities are also to be obtained. However, only the tangential 
velocities are considered unknown at this time. The evaluation of the 
meridional velocity is considered later, since its evaluation involves 
the solution of the flow equations in the meridional plane. In practice, 
both streamline and meridional plane equations have to be simultaneously 



satisfied, but to illustrate the analysis in the streamline direction, it 
is convenient to consider the solution of the streamline equations as part 
of an over-all iterative procedure. To a certain extent the developed an- 
alysis presented later integrates the two aspects of the solution of the 
flow field. Nevertheless, the selected over-all procedure still involves 
an iterative solution of the streamline positions* 

The analysis in the streamline direction is illustrated by con- 
sidering the design stations through a first stage. At inlet to the 
stage, station 0, the total pressure, P oa , the total temperature, ! O o > 
and the whirl angle, , will all be specified. The remaining quantities 
of interest at this plane will be determined from the values of meridio- 
nal velocity obtained from solution of the radial equilibrium and conti- 
nuity equations. The tangential component of velocity will be given by 


V Ro - v„ 0 

At station 1, the stator exit, the total pressure will be im- 
plicitly specified by the row total -pressure-loss coefficient; the total 
temperature will be equal to the streamline value at the stator inlet, 
and the tangential velocity will be explicitly specified or the flow 
angle will be given. 

The total -pressure-loss coefficient is defined as follows: 

Poo Po, (6) 

Pc, - P. 



Hence, 



( 7 ) 
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where 


( 8 ) 


fr - f/ 

Thus, the total pressure, Pot , can only be obtained when both and V»ti 
are known. Since will be dependent on Pet , a simultaneous solution 
for Poi and Vfcfj is required. 

In preliminary design investigations the tangential velocity 
will normally be specified, but as a design is finalized, it is conveni- 
ent to specify the flow angle. In this manner a particular variation of 
flow angles is obtained directly rather than by an iterative procedure. 
When the flow angle is specified, the tangential velocity is obtained 
from the following: 


14 '/ 

^ 9 c 


i4n*/ 7 

o J dp / O i J 


JL 


V^L/ V yyf I (isCS /\\ (9) 

If \flAi is specified, Equation 9 can, of course, be used to obtain the 
flow angle. 

Thus making the assumption that v* t wi 1 1 be determined from 
radial equilibrium and continuity considerations, the stator exit and ro- 
tor inlet velocity triangles can be completely defined using standard 
formulas and techniques. 

The next station, 2, is the stage exit. Here again the abso- 
lute values of total temperature, total pressure, velocity, and flow angle 
are considered the principal unknowns. None of these quantities will be 
directly specified. The total temperature will be obtained from the speci- 
fied work which is readily expressed in terms of the total temperature 
drop. Thus, 




( 10 ) 
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where W j s the work extracted along the streamline in Btu per Ibm. 

The tangential velocity component is obtained using the Euler 
work equation; that is, 

U, Vu.1 - 3-0 & T ° 


V, 


tta. 




( 11 ) 


where u, and are the blade speeds at rotor inlet and exit, respectively. 

The stage exit total pressure will normally be obtained from the 
rotor row tota 1 -pressure-loss coefficient,^, but for preliminary analy- 
ses it is convenient to specify stage isentropic efficiency, Is or rotor 
row i sent ropi c efficiency, , These last two opt i ons will avoi d the i t- 
eration which is involved when the rotor loss coefficient is obtained 
from a correlation within the computer program. For hand calculations 
of the velocity triangles, their use considerably reduces the time and 
effort required to complete a design. Stage exit total pressures are 
readily obtained from the standard definitions of tota 1 -to-tota 1 isen- 
tropic efficiencies. Thus, 

/ o / To A 


p« - p~ 


oo 


I - 
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OO 


JL 

X'l 


( 12 ) 


or 
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= Po. 
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1 o i T~ 
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(13) 

iressure-1 

oss coeff 

: icient is defined as 


Vr = 

Pox* 

'p'o-x. 

- pL 

- Pz. 


(14) 
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where r o 45 is the isentropic value of rotor exit relative total pressure. 
If there is no change in radius of the streamline from rotor inlet to 
exit, P OAS will equal P ot . However, if there is a change of radius, the 
isentropic relative total pressure is given by the expression 

x y 



+ - u. 
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Similarly, there is a relationship between the relative exit total pres- 
sure and the absolute pressure at stage exit which is readi 1y derived from 
the velocity triangles using the total -to-stati c temperature ratios. 


Pol _ t °Vt. 






To VT± 


z_ 

x-l 
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Z T To* 



(16) 


Thus, if Equation 14 is reexpressed as 


7, 


A ' 


P »*$!/ Pot \\ P<s*\ _ I 

1 ~(£J(W 


(17) 


fj/rz 

and the stage exit stat i c-to-total pressure ratio, ' , is expressed 

in terms of ar| d an explicit expression for /©* can be ob- 

tained from Equations 15, 16, and 17 in which i s the only unknown. 
Hence, as for the stator exit total pressure, the total pressure will 
have to be obtained from a simultaneous solution of P^x and 

When Po^y 7e>* , and have been obtained, the stage exit 

and rotor relative exit velocity triangles are readily obtained using 
standard turbine design techniques. 

In summary, the solution of the streamline flow conditions in- 
volve the tangential momentum or work equation and the pressure loss or 
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energy equation. In the developed form, the analysis uses the differen- 
tial form of these two basic equations as subsidiary equations in a simul- 
taneous solution of the radial equilibrium equation. The solution of 
these three equations provides a radial distribution of meridional ve- 
locity which satisfies the radial equilibrium requirement and the design 
specifications. For each station, the meridional velocity distribution 
is obtained for an assumed value of meridional velocity at the mean stream- 
line. The meridional velocity distribution, however, must also simultane- 
ously satisfy the mass flow continuity conditions. 

The Continuity Equation 

The location of the streamlines at each design station and the 
meridional velocity at each of the streamlines are determined from the 
mass flow continuity equation. For an axisymmetric flow passing through 
an axial section of an arbitrary annulus, the continuity equation is as 
follows: 

= 3-rr (, 8 ) 

'Sl 

where MX/- i s the total flow (lbm per sec), ^ is the static density (lbm per 
cu ft), and ^ and are the hub and casing radii of the station (ft). 

Throughout the analysis, total pressure, total temperature, and 
the tangential and meridional components of the absolute velocity are con- 
sidered as the principal variables. Hence, it is convenient to reformu- 
late Equation 18 in terms of these variables. Substituting for the den- 
sity using the standard formula, 


i 



14 


Equation 18 can be reexpressed as follows: 
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Since the radius dependent variables of Equation 20 in general will not 
be simple analytical functions of radius, the distribution of meridional 
velocity which satisfies continuity will have to be obtained using an it- 
erative numerical procedure. 

When the meridional distribution to satisfy the continuity re- 
quirement for the complete annulus flow has been determined, the location 
of the streamlines which satisfies the continuity requirements of the in- 
dividual stream filaments must also be computed. For the analysis, stream 
lines are selected so that any adjacent pair define an annulus containing 
a preselected constant fraction of the total flow. Thus, it is conveni- 
ent to introduce a mass flow function /£^A)which varies from zero to unity 
between the boundary streamlines. This mass flow function is defined as 
follows: , 
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J \/tyt C c$A 
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Hence, if the first streamline is at the hub and there are/ft streamlines, 

' J - 1 

the ^ streamline will have a mass flow function value of T~~T ' Thus, 

the radius of this st reaml i ne, Aj , will be obtained from the solution of 

the following equation 

J~t Po f f _ vj' + VjK 

~ X r R J To L Jgo'JCpT 

Thus, streamline positions throughout the flow field can be es- 
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S~t 


fj (22) 


tablished once the distribution of meridional velocity has been determined 



This distribution must simultaneously satisfy the continuity and radial 
equilibrium conditions for the design specification. 

The Radial Equilibrium Equation 

A fundamental assumption made in the analysis is that the flow 
is ax i symmet r i c . Hence, by equating the radial forces acting on a point 
P in the flow (see Fig 2), the condition for radial equilibrium is readily 
shown to be 

3l.4? = , 

^ cLA. si. si-bL. <sm (23) 

where Vj\, and are the tangential and meridional components of stream- 

line curvature, respectively. 

The third term on the right-hand side of Equation 23 results 
from a change in momentum in the meridional direction. However, in gen- 
eral, the change in both meridional velocity in the interblade row space 
and the sine of the streamline slope will be small. Hence, for the tur- 
bine design analysis, this third term will be assumed to be negligibly 
small. If it were decided that this term should be included in the analy- 
sis, it would be necessary to extend the analysis to consider the complete 
flow field defined by the trailing and leading edges of the upstream and 
downstream blade rows and the annulus walls rather than a section of the 
flow field on an axial plane. Thus, it would be necessary to make the 
assumption that the flow was axisymmetric at blade leading and trailing 
edges, which is considerably less valid than that concerning the axial 
symmetry at an interblade row design station. 

Omitting the third term on the right-hand side. Equation 23 
simplifies to 
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(24) 


This is the fundamental radial equilibrium equation used in the current 
analysis. The use of Equation 24 rather than Equation 23 simplifies the 
problem of solution considerably. However, the streamline slope and cur- 
vature in the mer i d i ona 1 p 1 ane , A and » still necess i ta te cons i de ra- 
tion of derivatives with respect to the axial direction, X. That is, 


-f 


and 


A - 

//, = - „ 


(25) 


( 26 ) 


Hence, the radiaT equilibrium equation contains derivatives with respect 
to both A and X. If the analysis were to be performed for an axi symmetric 
flow in an arbitrary duct, it could readi ly be extended to consider the 
path of individual streamlines in the meridional plane. However, in a 
turbine design-point analysis, it is unrealistic to assume that the axi- 
symmetric form of the radial equilibrium equation can be extended beyond 
the interblade row space into the blade rows. Thus, the boundary condi- 
tions for the meridional streamlines in the interblade row space are inde- 
terminate at the trailing edge and leading edge planes defining this space. 
Only the boundary streamlines, at the inner and outer annulus walls, are 
rigorously defined — by the assumption that these streamlines follow the 
contours of the annulus wall. In the absence of a rigorous analytical 
treatment for the slope and curvature of the flow in the meridional plane, 
it becomes necessary to adopt an arbitrary solution to the problem. In 
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the current analysis, streamline slope and curvature in the meridional 

plane will be treated in one of two ways: 

a. The slope and curvature of the boundary streamlines will 
be obtained from a definition of the wall contours, and 
then both A and VAm will be assumed to be linear functions 
of radius determined from values at the walls. 

b. The slope and curvature will both be specified arbitrarily 
as a function of radius. 

It is appreciated that other arbitrary solutions are possible. For exam- 
ple, existing design and off-design analyses developed at NREC and else- 
where for axial and centrifugal compressors have evaluated the slope and 
curvature of intermediate streamlines by spline-fitting curves, represent- 
ing these streamlines, through points at the design stations. For turbine 
design, however, it would be difficult to justify the increased complexity 
of the analysis, particularly when the influence of the blading on the 
streamline flow through the row is not considered. The radial and axial 
distribution Of loading and blade blockage throughout the rows will un- 
doubtedly have a larger influence on the values of /) and hi at a se- 
lected design station than the distribution of the flow at adjacent de- 
sign stations. 

With the incorporation of the assumption that slope and curva- 
ture are functions of radius, the radial equilibrium equation is of the 
form, 

\ ' 7Z = 7T - i f/t) < 27) 

Hence, the radial equilibrium equation can be solved station by station 
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throughout a turbine without an iterative procedure which involves other 
stations. The solution will depend only on the specified design require- 
ments, the preceding row losses, and the flow conditions at the preceding 
design station, and not on any data concerning the flow at subsequent de- 
sign stations. 

Since the continuity equation and the subsidiary streamline equa- 
tions have all been expressed in terms of the turbine design variables of 
total pressure, total temperature, and the velocity components, it is con- 
venient to express the radial equilibrium equation in terms of these vari- 
ables. Thus, it is necessary to reexpress the density and static pressure. 
The appropriate substitution for density has been presented earlier as 
Equation 19; the static pressure, P, is given by the standard relationship 



Hence, from Equation 19 and the differentiation of Equation 28, Equation 27 
can be reexpressed as follows: 
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Since the radial equilibrium equation has to be solved simul- 
taneously with two basic streamline equations (the tangent i a 1 -momentum 
equation and the tota 1 -pressure-loss equation), the principal analysis 
variables are considered to be > and -^7 

for its subsequent solution, Equation 29 is rewritten as 


. Therefore, 
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where C >f , (L u , 0 /3 , and C /v are coefficients which are readily obtained 
from the design specifications. It can be readily seen from Equation 29 
that C,| can be set equal to unity, equal to and that the remain- 

i ng coefficients will then be given by 


and 
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The development of the analysis, which consists essentially of establish- 
ing a calculational procedure for the solution of the flow field at each 
design plane is discussed later. I t wi 1 1 be shown that the two basic 
equations governing the flow in the streamline direction can be expressed 
as differential equations having the same form as Equation 30. 

I t wi 1 1 be noted that in deriving the particular form of the 
radial equilibrium equation to be used in the computer program solution 
of the flow field, it was assumed that the specific heat and the specific 
heat ratio, and X, are independent of radius. Similarly, throughout 
the discussion of the streamline equations it was implicitly assumed that 
these gas properties were known constants. For the analysis i t wi 1 1 be 
assumed that the specific heat will be constant at any axial station but 
may be specified to vary through the turbine. 
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MODIFICATIONS TO THE STREAM-FILAMENT APPROACH 
FOR VARYING SPECIFIC HEAT, MIXING, AND COOLANT FLOWS 


I nt roduct i on 

In the preceding section of this report in deriving the stream- 
line equations, it was assumed that the specific heat was constant. How- 
ever, the analysis is to form the basis for a computer program to be used 
for multistage and multispool turbine configurations. Thus, in many design 
investigations the variation in temperature through the machine will be 
sufficiently large that the variation of specific heat will be a signifi- 
cant factor. The analysis will be modified to accommodate a station-to- 
station variation of specific heat, and these modifications are presented 
in this chapter. 

The stream-filament approach assumes the flow field to consist 
of concentric surfaces of revolutions defined by selected streamlines. 

In an actual turbine it is unrealistic to assume that there will not be 
any mass transfer between these stream filaments. While the flow, which 
remains in the free stream throughout its passage through a turbine, will 
probably remain within its original filament, the flow which is affected 
by the viscous forces near blade surfaces will tend to migrate to other 
filaments under the influence of the complex static pressure field which 
is set up within a turbine. Experimental test data from turbine stages 
suggest that for many turbines this mixing is sufficiently severe that 
the basic stream-filament approach should be capable of modification to 
represent, at least qualitatively, this mixing. 

An additional factor which must be considered in many turbine 
applications Is the addition of coolant flows to the main stream. Coolant 
flows can be added to the main stream in a variety of ways. For example, 
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disk and rotor root coolant flow will enter along the hub contour, rotor 
convective cooling flows will normally enter the main stream at the outer 
casing, and where film or transpiration cooling is used, the coolant flow 
will be distributed across the annulus. Irrespective of the manner of 
admission, the principal effects which should be considered in any modifi- 
cation of the fundamental stream-filament approach will be to increase the 
flow and to change the temperature level of the main stream. The modifi- 
cations which are made to accommodate a specified coolant schedule in the 
analysis are discussed, in the last section of this chapter. 

Variations of Specific Heat 

The simultaneous solution of the radial equilibrium and conti- 
nuity equations at the turbine design stations is considerably simplified 
by the assumption of a radially constant value of specific heat. Using 
the standard formulas for tota 1 -to-s ta t 1 c temperature and pressure ratios, 
which assume a constant specific heat, the form of these two basic equa- 
tions are amenable to solution in terms of the principal variables /?, To , 
V/ k , and Thus, the principal justification for the use of radially con- 

stant specific heat and specific heat ratio must be that it considerably 
simplifies the numerical solution of the radial equilibrium and continuity 
equations. However, it should be noted that at a stator exit station of 
a design In which there is a zero reaction at the stator casing section 
and zero reaction at the following rotor hub section, the radial variation 
of specific heat will be as large as that which exists across a stage. 

When two or more design stations of differing specific heat are 
used in the foundation of the streamline equations, an appropriate mean 
value of specific heat will be used. The work, blade row loss, and 
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i sentropic efficiency expressions all involve conditions at two stations. 
When the specific heat is not assumed to be constant, the Euler work equa- 


tion has to be expressed as 
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3. r 

) 

S*P dT o ~ ^ ~ ) 


(33) 


(34) 


where 


g p s 


C^I + c 






Similarly when the equations involve an isentropic expansion from one sta- 
tion to another, the expansion index will be based on a mean value of spe- 
cific heat. Thus, for a stage pressure ratio, the isentropic 

temperature ratio will be defined as 
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The approach of using an appropriate mean value of the specific heat when 
the thermodynamic relationship relates state conditions at differing sta- 
tions will lead to more complex expressions for all the streamline equa- 
tions. For example, the stage efficiency equation wilt contain three 
specified values of specific heat. That is, 
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For a single stage the actual variation between C^, C^, , and will be 
relatively small, and the calculated efficiency will differ little from 
a value obtained from a design in which a constant specific heat has been 
assumed. In a multistage unit, the over-all variation in may be more 
significant. Thus, it could be argued that the ideal work should be based 
on individual stage isentropic temperature drops calculated from the ap- 
propriate stage values of pressure ratio, stage mean specific heat, and 
inlet total temperature. However, the isentropic temperature drop for 
the complete spool, which will be used in the expression for over-all ef- 
ficiency, will be simply obtained from the over-all pressure ratio, the 
spool inlet total temperature, and a mean specific heat ratio for the 
spool , 


Mixing 

In the current analysis, no attempt is made to accurately repre- 
sent the complex viscous flow phenomena which occur in a turbine, nor to 
assign any loss in blade row efficiency to the interfilament mixing pro- 
cess. Experimental investigations using flow visualization techniques and 
the detailed analysis of turbine performance using data from radial and 
circumferential surveys clearly show that mixing can be sufficiently se- 
vere that the fundamental streamline analysis should be modified to repre- 
sent at least qualitatively the mass flow mixing which occurs between the 
axisymmetric stream filaments. 

While secondary flows have associated losses, and one of the mix- 
ing mechanisms is undoubtedly these secondary flows, the correlation of 
losses to be used in the design analysis program has loss levels which 
are in excess of those due to profile loss. Hence, for the purpose of 
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this analysis, no loss is assumed to be directly associated with mixing. 
Thus, if a total pressure or total temperature profile is modified, the 
new profiles will have the original mass flow weighted values. This con- 
sideration leads to a formulation of a mixing model of the following form: 


PoJ = i- Xj l x * (37) 




and 
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where A oj and /Cj are the mixed values of total pressure and total tempera- 
ture for streamline^. The mixing pa rameter ,)£j » can be specified for indi- 
vidual streamlines, but the same value will be used for both total pressure 
and total temperature. Since the TfV stream! i nes will def i ne(7i-j) st ream fila- 
ments of equal flow, the flow associated with individual st reaml i nes , , 

will equal for internal streamlines and for the hub and casing 

boundary streamlines. Both equations satisfy the set requirement that the 
mass flow weighted mean values of mixed and unmixed profiles will be equal. 
With the above formulation for mixing, any radial distribution of the stream- 
line mixing factor may be specified with the exception of those for which 

No mixing Is, of course, one of these special cases, but then it 
will be unnecessary to use Equations 37 and 38 since the streamline values 
of total pressure and temperature will be unmodified. If experimental data 
Indicate that complete mixing occurs near the annulus walls but that virtu- 
ally no mixing occurs in the central portion of the annulus, the appropri- 
ate streamlines could be assigned values of 1.0 and 0.0, respectively. 


Even though the mixing parameter has been selected as a stream- 
line dependent variable^ it can, of course, be specified as a constant for 
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a particular row. In these circumstances, Equations 37 and 38 simplify 
to become 

Poj =//-xJP oj + X-K (39) 

and ^ ’ (I-X) loj + X-T; (40, 

where /*> and 7^> are mass flow weighted values for the original profiles. 

In many analyses it is probable that constant values of X wi 1 1 
be specified. For a fully mixed solution with X— U° at each streamline 
the mixed values would equal the corresponding mass flow weighted value. 

Having chosen a mathematical formulation for the mixing, it is 
necessary to decide at what point it is to be introduced into the analy- 
sis. Since mixing is related to the flow within a blade row, it would 
appear logical to specify mixing parameters for blade rows rather than 
design stations. Thus, the mixing will occur before or after the design 
stations at which radial equilibrium and continuity equations are satis- 
fied. The modification from streamline to mixed values of absolute total 
pressure and total temperature will be made downstream of the plane at 
which radial equilibrium and continuity equations are satisfied and will 
be used as inlet values for the following blade row in which mixing is 
assumed to occur. Thus, revised values of inlet total pressures will be 
used in the total -pressure-1 oss coefficient expression, and the revised 
values of inlet total temperature will be used as the datum for the tem- 
perature drop through the rotor row. The selections of the row inlet 
rather than the row exit plane considerably simplifies the numerical cal- 
culations at the following interrow design station where radial equilibri 


and continuity are again satisfied. 

Since the mass flow weighted values of the revised total pres- 
sures and total temperatures will be unchanged from the corresponding 
stream-filament values, the stage mean values of efficiency and work out- 
put will not be directly affected by specified mixing. Some small effects 
on mass flow weighted stage efficiency is to be expected, since the levels 
of loss coefficient will be dependent on the computed blade geometries. The 
distribution of flow angles and velocities will, of course, depend on the 
solution of the radial equilibrium and continuity equations at the follow- 
ing design station, and hence the modification of the radial gradients 
of enthalpy and entropy due to mixing will have some effect on the levels 
of loss. Unfortunately, i t wi 1 1 be necessary to develop the computer 
program before the influence of mixing can be fully investigated. Nev- 
ertheless, it is clear that significant effects are to be expected if 
large gradients of total pressure and total temperature (which result 
from the computed radial distributions of loss coefficients and speci- 
fied work extraction) are assumed to fully mix. Since some of the total 
pressures used in the tota 1 -pressure- loss equation will be greater than 
their corresponding unmixed values, it is quite possible that stream- 
filament efficiencies, based on mixed values of total pressure and total 
temperature will be in excess of 100 per cent in some instances. Such 
anomalies have been observed experimentally in turbine stages tested with 
a radial distribution of stage inlet total temperature. In order to pro- 
vide guidance to the program user on the effects of the specified mixing, 
both the streamline and mixed values of total pressure and total tempera- 
ture will be included in the computer program output. It should be em- 
phasized however that the particular values of the mixing parameter which 


27 



will be specified for a selected turbine design-point analysis will have 
to be fully investigated using the computer program as an analytical tool 
In conjunction with detailed experimental data for similar turbines at a 
later date. 

Coolant Flows 

For some high inlet temperature turbines, the amount of coolant 
mass flow admitted to the main stream will be sufficiently large that sig- 
nificant mismatching of stages would result if the coolant mass flow sched- 
ule were not an integral part of the design-point analysis. The tempera- 
ture level of the coolant must also be considered. The relevant tempera- 
ture is, of course, the initial temperature of the coolant which will often 
be that of the compressor delivery air. Even though at the point of admis- 
sion the coolant flow, which has been used for blade cooling may be close 
to the temperature of the main stream, it will have attained this tempera- 
ture because heat has been transferred from the main stream by way of the 
% 

blading. Hence, the local temperatures to be used in the design-point 
analysis should be derived from consideration of the heat balance equa- 
tion. Thus, coolant mass flows and temperatures will have to be speci- 
fied for the design-point analysis of a cooled turbine. However, since 
the coolant may be admitted in a number of ways, it is inappropriate to 
consider the pressure level of the coolant, its flow direction, or its 
radial distribution in a general design-point analysis. Similarly, any 
losses associated with the introduction of coolant cannot readily be ac- 
commodated within a stream-filament analysis other than by the program 
user directly specifying row loss coefficients or additional loss factors 
for the rows whose performance will be affected by the coolant. 
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As with the mixing previously discussed, the addition of cool- 
ant will be considered as a blade row rather than an interblade row pro- 
cess. Hence, the coolant flow schedule will be specified in terms of 
the mass flow and temperature of the coolant for each of the cooled rows. 
The treatment of the two aspects of the coolant flow is discussed, in 
turn. 

Considering the coolant mass flow which is added to the main 
stream in the blade row which precedes design station/H, the continuity 
equation for this station will be written as 

Ac 

^ * X, * Ajur = IT j *■ <**■ m 

*-k. 

where A^is the coolant flow added to the main stream in the blade row 
preceding station syi* 

In the calculation of the gas temperatures resulting from the 
addition of coolant to the main stream, it will be assumed that the cool- 


ant will be uniformly distributed between the main stream filaments. 

Hence, the temperatures will be modified on a streamline basis with the 

th 

new temperature, / oj , of the j streamline being given by the expres- 
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With the addition of cooling being considered a blade row phe- 
nomenon, the change in temperature is assumed to occur between design 
stations. As with the mixing previously discussed, a decision has to be 
made whether to make the temperature level correct at a blade leading 
edge plane (which is immediately downstream of a design axial station) or 
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at a trailing edge plane ( i mmedi atel y upstream of the next design sta- 
tion at which radial equilibrium and continuity are to be satisfied). 

For stator blading, where streamline total temperatures are constant 
through the row, the final temperatures would be identical for the two 
alternative methods. For the rotor, across which there is a change in 
absolute total temperature level, there will be small differences be- 
tween the results of the two methods because the radial gradient of tem- 
perature will be changed in any design where radial variations of work 
output has been specified. However, the numerical solution of the de- 
sign problem and construction of the computer program are simplified by 
correcting the temperature level at the inlet plane of a cooled row; 
this approach has been selected. 

For a stage in which coolant is added to the main stream in 
either or both of the blade rows, there is no generally accepted defini- 
tion of the stage isentropic efficiency; the expansion process is, in 
fact, no longer isentropic even in the absence of losses. For the cur- 
rent analysis it is assumed that all the flow which leaves the stage 
contributes to the work output of that stage. Thus, if station/t^is a 
stage exit, the stage total temperature drop,^*o, along any selected 
streamline is obtained from the specified total power output (in Btu’s) 
and distribution of output and the total flow,/^T. That is, 
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The total temperature "ToJ is obtained by correcting the stage inlet 
total temperature for both the stator and rotor coolant f 1 ows , 
where 
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(44) 




I t wi 1 1 be assumed that the ideal power output of the stage is 
that available from the expansion of both the main flow and the coolant 
flows across the stage pressure ratio. That is, 


From Equations 43 and 45 the expression for stage efficiency becomes 


(45) 



where 
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The over-all efficiency of a spool having cooled stages will be expressed 
as a summation of actual power outputs divided by a summation of ideal 
power outputs. The summation of actual power outputs is, of course, 
directly obtained from the design specifications; the ideal power out- 
put of the turbine inlet flow and that of the individual stage coolant 
flows will be evaluated using the appropriate values of pressure ratio 
and inlet total temperature. 


Concluding Remarks 

Relatively simple modifications to the constant specific heat 
stream-filament analysis are used to take account of variations of spe- 
cific heat, interfilament mixing, and the addition of coolant flows. 
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There are two principal justifications for the use of relatively simple 
models. Firstly, they are readily accommodated within the structure of 
a purely s tream-f i lament analysis. Secondly, for the majority of appli- 
cations the effects will be relatively unimportant. For many single or 
two stage turbines the actual variation of specific heat will be almost 
insignificant. Although the effects of interfilament mixing cannot, at 
this time, be fully investigated, it would appear probable that only for 
designs in which large gradients of work extraction are specified will 
the total pressure and total temperature profiles be greatly changed by 
mixing effects; however, experience with stream-filament designs indi- 
cate that if too severe a radial gradient of work extraction is speci- 
fied, the blading geometries become mechanically unacceptable or it be- 
comes impossible to satisfy the requirement of radial equilibrium. Total 
coolant flows will, in general, be a relatively small percentage of the 
machine inlet flow and in multistage applications will be limited to the 
early blade rows; while the model for the addition of coolant to a par- 
ticular stage is relatively simple and open to question, the correction 
to the levels of flow and temperature for the calculation of the geometry 
of the downstream stages will be valid. The avoidance of stage mismatch- 
ing is the most important aspect of cooled turbine designs, and this is 
adequately considered in the proposed model. 

While mixing will invalidate the concept of a stream-filament 
efficiency, the mass flow weighted stage efficiencies will not be di- 
rectly affected in the mixing. All efficiencies will be influenced by 
the introduction of coolant. While the selected definitions of effi- 
ciency for cooled stages are not necessarily generally accepted, the 
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program output will contain all the total pressures and temperatures gen- 
erated in the design-point analysis; from these data, alternatively de- 
fined efficiencies may be readily obtained. 
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DEVELOPMENT OF THE LOSS CORRELATION 


I ntroduct i on 

It is self-evident that the most sophisticated numerical analy- 
sis of a selected turbine design requirement would be of little value if 
blade row performance data used in conjunction with the analysis were not 
consistent with the blading to be used for the design. In adopting a 
stream-filament approach to turbine design, the loss data required are 
for elements of blading in a stage environment and as such are not readily 
available. Simple cascade data can provide a guide to the expected per- 
formance, but profile loss is often a relatively small percentage of the 
total loss. Apart from the basic question of loss level, cascade data 
are frequently presented in terms of the over-all profile geometry such 
as blade angles, solidity, and thickness/chord ratio but are rarely re- 
lated to the detailed design of the profile; the latter is undoubtedly a 
significant factor, since it ultimately governs the growth of boundary 
layers. For a design-point analysis, it is desirable to relate loss 
levels to the design requirements which can be expressed in terms of the 
over-all loading of the section required by the design. A basic assump- 
tion made in the present analysis is that the blading is suitably designed 
to avoid the excessive losses associated with unsuitable profiles, solidi- 
ties, or incidence settings. Since every turbine is a compromise between 
aerodynamic and mechanical requirements, it is inevitable that blade row 
losses are often increased as a result of mechanical requirements Hence, 
the datum level of loss which is to be an integral part of the computer 
program must in some instances be factored to compensate for any antici- 
pated increase in loss. 
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Since there is little a va i 1 a b 1 e da ta on wh i c h to base e 1 emen t 
losses, the present loss correlation was obtained from a correlation of 
achievable stage efficiencies. The manner in which the correlation was 
obtained and some checks on the validity of the correlation are presented 
in the following sections. Loss factors and kinetic-energy-loss coeffi- 
cients are also briefly discussed. 

Data for the Loss Correlation 

Since the loss coefficient data are to be applied along stream- 
lines passing through blade rows, it would appear that the correlation 
should be based on experimental traverse data from turbine stage tests. 
Unfortunately, very few detailed radial traverse investigations are under- 
taken during turbine stage performance testing and even less data are 
readily available in the published literature. Hence, it is necessary 
to base the correlation on over-all stage results using the frequently 
adopted mean-line approach to turbine analysis. While considerably more 
data exist for complete stage designs, it is important that these data 
reflect what is achievable in terms of efficiency level rather than that 
which is measured when a design is inefficiently bladed or is operating 
at off-design conditions. 

The principal source of data used for the correlation is the 
stage efficiency correlation of Reference l« This reference contains a 
contour plot of achievable stage efficiencies correlated against the 
over-all stage design parameters of stage loading factor , ATo /tl^ 

and stage flow factor, This particular plot, which is based on 

the results from a large number of turbine stage performance tests, is 
reproduced as Figure 4. The efficiency levels shown are corrected to 
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remove that loss which is believed to be due to tip clearance effects. 
l n general, the shape of the contours and levels of efficiency shown in 
the correlation of Reference 1 differs little from similar correlations 
obtained by other ai rcraf t gas turbine manufacturers. The stage load- 
ing factor and stage flow factor do not uniquely define the turbine de- 
sign velocity triangles but by making four assumptions, the stage perfor 
mance correlation can be used as a source of individual row tota 1 -pressure- 

loss coefficients. 

These assumptions are: 

1. Rotor and stator loss coefficients will be equal when their 
relative design requirements are identical. 

2. The turbines presented in the correlation were designed 
with 50 per cent stage reaction at the mean line. 

3. The axial velocity is constant through the stage. 

4. The stator exit Mach numbers were 0.8. 

Each of these assumptions is considered in turn. 

* 

The first assumption is a fundamental requirement of a stream- 
filament analysis; one of the basic requirements of the correlation is 
that it can be used for any design section irrespective of whether it is 
a rotor or stator row and its radial location in the annulus. In the 
event that use of the analysis program, in conjunction with test results 
from turbine stages designed using the stream-filament method, provides 
basic data which do not support this assumption, loss correction factors 
will have to be employed to adjust loss levels. 

The second assumption is expected to be valid. The majority of 
turbines used in the correlation will have mean section stage reactions 
(defined as the ratio of static temperature drop across the rotor to the 
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stage total temperature drop expressed as a percentage) in the range from 
40 to 60 per cent. Even though the stages may not in every case have a 
50 per cent stage reaction, the effect on efficiency level of any small 
difference from 50 per cent is likely to be small. Although not a neces- 
sary assumption for the derivation of row loss data, it is of interest 
to note that the majority of turbines used in the correlation would have 
had rotor hub-section reactions sufficiently far from impulse conditions 
that the stage efficiency levels would not be greatly affected. 

From published engine information, it is evident that the tur- 
bine stages used in the correlation were designed with significant amounts 
of annulus flare, and from this it can be concluded that many of the de- 
signs will closely approximate constant axial velocity designs. Hence, 
the third assumption is reasonable. 

The actual Mach number levels for the stages tested are unknown, 
but it is extremely likely that the assumption of a stator exit value of 
0.8 is reasonably valid; deriving stage pressure ratios for the selected 
velocity triangles shows these lie in the range 1.5 to 2.0:1, which is 
the range of stage pressure ratio likely to be used in aircraft engine 
appl i cat i ons. 

With these four assumptions, complete velocity diagrams may be 
obtained for any turbine design point which is defined in terms of its 

stage loading factor, T~ a , and stage flow factor, S ~ . 

Hence, using efficiencies from Figure 4^ it is possible to deduce the row 
total pressure loss coefficients of selected turbine designs. From the 
definition of reaction 

R = JL zZk = / * -2 - V* < 47 > 

T 0 „ -7^, w 
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the ratio ’ K, /M is obtained. Hence, the stator exit angle^is readily 
computed from the expression 


(48) 


(49) 


The relative flow angle, , is given by the equation 

//= 

With the angles of the stator exit and rotor relative inlet diagrams de- 
fined by selected values of ^and <j>, these diagrams can be completed in 

terms of the velocity parameters/-^ Jj-=. etc.) by assigning a value 

J VToi t ' 

to the stator exit Mach number. 

The assumptions of 50 per cent reaction and constant axial ve- 
locity across the rotor produce stage exit and rotor relative exit dia- 
grams which are similar to those at stator exit; that is, and 

^ - -fy . To completely define these stage exit velocity triangles, 
the ratio of stage exit to inlet total temperature, l(,x Alo’ is rec l uired 
and this ratio is simply obtained from the definition of stage loading 

factor. That is, 

To, _ . -y. 
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Thus, all the relevant velocity parameters can be computed using the blade 
speed parameter at stage exit ^T CA _ where 17 , = and 

U A 

In addition to the normally computed velocity triangle quanti- 
ties, the stator exit absolute and rotor exit relative values of total -to- 
static pressure ratios, P °/p, and f? , are also required in order 
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to derive row total pressure loss coefficients. The relationship between 
stage isentropic efficiency,^, and row loss coefficients, X/ anc * ^ » 
i s as fol lows: 



The above expression is an approximation, but it is of a high 
order of accuracy when the efficiency level is in excess of 85 per cent. 
Using the appropriate values of efficiency obtained from Figure 4 for the 
selected design points, Equation 52 can be solved for the rotor loss co- 
efficient,^, provided the stator loss coefficient,^, is also known. 
Values of Y/v are obtained from design points having a unity stage loading 
factor; for these turbines the rotor inlet angle is zero (see Equation 49) 
and the two row loss coefficients are assumed to be equal, making it pos- 
sible to evaluate them from Equation 52. 

The loss data used for the correlation are based on four stator 
configurations, these having exit angles of 71 * 8 , 60, 51*3, and 45 de- 
grees. Additional data points were obtained for various rotor blades by 
considering design points having these stator exit angles. These design 
points lie on four straight lines In the diagram. The equations of 

these lines are obtained from Equation 48, which can be reexpressed as 

±. - }<^ <f> - 3 - 5 “ (53) 

where the constant K of any selected line is equal to the tangent of the 
stator exit angle, The loss coefficient correlation is based on 

18 design points that span the area of the efficiency correlation which 
is substantiated by measured turbine performance. (Turbine test data 
points have been omitted from Figure 4 but are shown in the corresponding 
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figure of Reference 1.) 


Correlation of Total -Pressure-Loss Coefficients 

One of the widely used correlations (that due to Soderberg) 
assumes loss is principally a function of deflection. Hence, one of the 
first correlations attempted was tota 1 -pressure- 1 oss coefficient versus 
row deflection shown in Figure 5* Points having the same exit angle are 
joined by lines in this figure, and values of row reaction are noted 
against individual points. It will be seen that loss coefficients gen- 
erally increase with deflection but that row reaction has a strong influ- 
ence on the level. The same data points are replotted in Figure 6 using 
row reaction as the abscissa. (Row reaction is defined as I ~~ Vpy^) 

where and Vi^, are the inlet and exit velocities relative to the blade 

section.) In Figure 6 the values of deflection are noted against the 

C 

points. Assuming row def 1 ec t i on , u , and row velocity ratio, , are 

the relevant parameters, the data are correlated to a reasonable degree 
of accuracy by the following expression: 

y - (0,cco7^ i-0-COOOO5'& )^ 3 ' ? lo ( f 

where <5" is the row deflection in degrees. 

While the above correlation reproduces the efficiency contours 
of Figure 4 with acceptable accuracy, it cannot be considered acceptable 
from theoretical considerations. The above correlation implies that when 
the velocity ratio, , falls below 0.45, the level of loss for a 

selected deflection begins to increase. This result is unacceptable, since 
it is to be expected that losses will decrease smoothly as the over-all row 
acceleration increases at a constant deflection. The high reaction blading, 
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wi th low values of 


/ 14*. will have high values of exit angle. Hence, 
the possibility that the increase in loss was due to trailing edge block- 
age effects was investigated. In Figure 7 the loss coefficient data were 
replotted with blade row exit angle as the abscissa. Even though the in- 
dications are that a theoretically acceptable variation of loss with re- 
action could be obtained if the exit angle were introduced as a parameter 
in the correlation, the over-all correlation would be complex. 

In the initial correlation, the inlet and exit angles appear 
together as deflection; a review of the basic data showed that the unac- 
ceptable form of the reaction factor could be traced to a failure to dis- 
tinguish between blades having the same deflection but significantly dif- 
ferent design requ i rements. For example, in Figure 6 a 70 degree deflec- 
tion can be located on the three curves having exit angles of 51.3, 60, 
and 71.8 degrees. The inlet angles corresponding to 70 degree deflec- 
tions are 18.7, 10, and -1.8 degrees, and the row reactions are 0.34, 

0,49, and 0.69, respectively. It will be seen that the loss coefficients 
of the first arid last are greater than that of the intermediate blade. 

The use of deflection as a correlating parameter is clearly the reason 
for the theoretically unsound dependence on velocity ratio given in Equa- 
tion 54. The basic difference between the three profiles considered in 
the above example is the change in angular momentum that is required of 
the profiles. The tangential loading of the profile is directly related 
to the change in tangential velocity across the row, and hence it is logi- 
cal to consider tangents of the inlet and exit angles rather than the 
actual angles In the correlation. 

If the inlet and exit axial velocities are assumed to be equal, 


the change in tangential velocity is proportional to - +• Uv'/W 

for stators and Kvvy\„ - for rotors. Thus,| b.n 1 can 

be selected as a tangential loading function. To investigate a possible 
correlation of the form 

Y= (55) 

loss coefficients were divided by their tangent function and plotted 
against the row velocity ratio (see Fig 8). It will be seen that the de- 
pendence of loss level on the velocity ratio is now of the expected form, 
with loss level increasing as the design requirement approaches impulse 
conditions, 1 0 . There would appear to be some dependence on 

exit angle which again indicates the possibility that trailing edge block- 
age might be a factor in the level of loss. Since blade trailing edge 
thicknesses are frequently selected to be a constant fraction of the blade 
pitch or chord, it might be expected that trailing edge loss would in- 
crease with trailing edge exit angle. The trailing edge blockage is rep- 
resented by the ratio of trailing edge thickness to the exit flow area, 
HvfWW) . Hence, if the ratio is approximately constant for 

standard blade design practice, the level of loss might be expected to be 
a function of . Based on the loss data, a correlation of the type 

Y = 1 C/V*\ (56) 

~ (c-u 7 v V«a) 

was investigated. The data for a correlation of this type are plotted in 
Fi gure 9- 

Considering Figures 8 and 9 detail, a number of points emerge. 
Firstly, both figures exhibit the increase of loss with velocity ratio 
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which might be expected (in contrast to a correlation using deflection). 
Secondly, there is some reduction in scatter when the correlation includes 
a correction to represent possible trailing edge blockage losses. Thirdly, 
the scatter shown by the points beyond a velocity ratio of 0.55 makes the 
selection of a final form of the loss correlation extremely difficult. 

In general, the data show not an unreasonable nor an entirely unexpected 
result. With a velocity ratio of 0.5 or less, it is possible to design 
blade sections with little or no additional loss d,ue to local surface ve- 
locity diffusions. As the velocity ratio increases towards its impulse 
value, it becomes increasingly difficult to avoid an increase in loss due 
to suction surface diffusion effects. The actual level of loss will depend 
on the actual severity of the adverse static pressure gradient, which will 
in turn depend on the detailed design of the section. In poorly designed 
blades, additional losses may also result from increased diffusion on the 
initial portion of the pressure surface. Thus, not only is an increase in 
loss level to be expected as the velocity ratio increases, but some scatter 
might be expected in the data used for the correlation as a result of 
varying standards of blade design. A theoretical study of row losses pre- 
sented In Reference 3 relates the total blade surface diffusion,!^ , to 
reaction (defined as I ) and tends to substantiate the increase 

in loss level for velocity ratio above 0.55. Figure 10 of Reference 3 
shows a zero value of total diffusion at a velocity ratio of 0.45 and an 
almost linear increase as the velocity ratio is increased to unity. 

The data used for the current loss correlation study is limited 
to a velocity ratio of 0.72. However, for a general application of the 
correlation, i t wi 1 1 be necessary to extend the correlation beyond the 
range of the basic data. Since the characteristic form of the data will 
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make an extrapolation difficult, a data point corresponding to a blade sec- 
tion beyond impulse has been added to those initially derived. This point 
was obtained from stage performance data of a high-pressure-rat io turbine 
presented in Reference 2. The point, which has been included in Figures 
8 and 9, is for a high deflection rotor operating beyond impulse. This 
data point tends to confirm a correlation of the type suggested in Equa- 
tion 55 rather than Equation 56. However, the data points corresponding 
to velocity ratios of approximately 0.71 originate from an area of the 
stage efficiency correlation where there are a few test results. These 
loss coefficients were derived from points on a line (^- 2 ^- 1 ) which 
passes through the efficiency contours of Figure 4 where they are most 
open to question. 

While a precise correlation is difficult to select, the over-all 
form of the correlation can be chosen with a fair degree of confidence. A 
good correlation for the function of velocity ratio would appear to be two 
straight lines; below a velocity ratio of 0.5 the level of loss is almost 
independent of row reaction but beyond 0.5 the level of loss increases 
approximately* 1 inear ly with velocity ratio. However, because the computer 
program will use a forward-stepping procedure involving the derivative of 
the loss coefficient, a correlation continuous in the first derivative of 
loss coefficient with respect to velocity ratio was selected. In addition, 
due to the lack of certainty in the selection of a correlation, the co- 
efficients of the correlation will be input quantities for the computer 
program. The selected correlation is of the form 



1 Un ~ /Ua| 



if 


\L 


Oa 
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and 


(57) 


Y , 1 la W\^ - j a?>/k>xl ( a t + a 7 M*') < f ^ 4 a , 

(<x r -t-Qi oca /i *.*) L '^W J V** 

where Q f to Qy are coefficients which may be selected by the program user. 
The veloci ty-ra t i o-dependent function using a t = 0.055, ft?. = 0.15, c* 3 = 

0.6, Q£ - 0.03, Qj = 0.157255, and #£ « 3.6 is shown in both Figure 8 and 
9. This particular type of correlation is suitable for a design analysis. 

If the coefficients a,, q 2 , c\^, and are selected initally coefficients 
o 7 and can be readily evaluated from the conditions that the values of 
Y and ^^(^o^are continuous at a value of 

Since the value of Y computed from Equation 57 can become ex- 
tremely large as an exit angle approaches 90 degrees, the correlation will 
also include a maximum value of the tota 1 -pressure-loss coefficient, 

Mean-Line Stage Performance Prediction Utilizing Loss 
Coefficient Correlation 

Although the form of the loss coefficient correlation is theo- 
retically acceptable, the values to be assigned to the constants are dif- 
ficult to determine with any degree of certainty. In order to arrive at 
the final correlation to be recommended for the computer program, a number 
of stage performance predictions were undertaken. The principal objective 
of this investigation was to establish the values of the constants that 
gave the best over-all agreement with the achievable efficiency carpet plot. 
Predictions were undertaken initially at four levels of stage loading fac- 
tor. The results of these efficiency predictions are shown in Figures 10, 
11, 12, and 13 for values of of 1.0, 1.5, 2.0, and 2.5, re- 

spectively. Each figure shows as a reference the experimental data corre- 
lation and predicted curves for three alternative forms of loss coefficient 
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correlation. The correlations compared are 
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( 58 a) 


Y=|t^-U / ^4^«+0; l |^-H + c ’ l l^-°^ <58k) 

A] (58c) 

these alternative correlations assume, for convenience, two linear branches 
of the velocity ratio function rather than the more complex form of Equa- 
tion 57 * (in effect, coefficient O-j is set equal to zero and the condition 
for continuity in the derivative of Y with respect to is ignored.) 

Each of the alternative correlations has regions in which it fits the ex- 
perimental data more closely than do the others. All produce a contour 
plot with the peak efficiency ridge following the 60 degree exit angle sta- 
tor design line. it is in fact very difficult to decide at this time 
which of the three is the better correlation. Efficiency contour plots 
were prepared using correlations of Equations 58 a and 58 b, and these are 
shown in Figures 14 and 15. These two correlations differ only in the 
correction for trailing edge exit angle. In Figure 14 the trailing edge 
blockage correction is included, and i t wi 1 1 be seen that, at any selected 
loading level, the efficiency falls more rapidly as the flow factor is de- 
creased below its optimum value than when it is increased. In Figure 15, 
in which the correlation of Equation 58 b is used, the contour plot becomes 
more symmetrical about the peak efficiency ridge. Unfortunately, if these 
figures are compared with the original data (see Fig 4), it would appear 
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that the correction should be applied for large exit angles but omitted for 
small angles. Some improvement could be obtained by changing the constants 
in the term but a review of data on which the correlation is 

based suggests that this is unwarranted at this time. Data points are 
shown in Figure 8 of Reference 1, and from this plot it can be appreciated 
that alternative positioning of the contours could be justified. An addi- 
tional factor to be considered is that the turbines used in the correlation 
cover a considerable time span, and thus it is almost inevitable that de- 
sign standards for both stage aerodynamics and blading have changed in 
this period of time. 

In summary, a correlation of the form given in Equation 57 will 
be suitable for a design analysis program. However, since there is some 
uncertainty in what are the appropriate values of the coefficients, the co- 
efficients a, to will be made part of the input specifications for the 
computer program. The recommended correlation is as follows; 

Y - 1. Aw _ I C-C55" *- C 'S'/ ^ - 0 i)\ i f Vw\ > Q . I 
(CL l G-tf uT> f\**) l 'V** ' ) 

and Y-= 1 kw/U - fc-pT, -I- k if &<C-b (59) 

/ , Y. / \ ' vw J VuX 

(C’b f-0 ^ 

However, this particular correlation is based on a mean-line analysis of 
turbine stages which would have efficiencies consistent with the n achiev- 
able M efficiency correlation of Reference 1. Hence, it is quite possible 
that it will be necessary to modify the correlation at some future date. 

For example, the values assigned to coefficients and of Equa- 

tion 57 may be found to be pessimistic for well-designed blades in which 
excessively high suction surface Mach numbers and large diffusion gradients 
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are avoided even at high values of velocity ratio • 

Stream-Filament Prediction of Stage Performance 

In arriving at the correlation, mean radius design analysis has 
been used with over-all stage efficiencies. However, the stage results are 
for blading which varies considerably with radius with respect to the cor- 
relating parameters. In addition, the correlation is to be used on an ele- 
ment basis in the design analysis. It is therefore necessary to investi- 
gate what effect the use of the correlation has on predicted stage effi- 
ciency when it is applied using the stream-filament approach. For this 
purpose two designs were selected and were analyzed using an existing NREC 
stream-filament design program. For the investigation tota 1 -pressure- loss 
coefficients were based on the correlation of Equation 58a. The designs 
were arbitrarily selected; one can be considered typical of the first tur- 
bine stage of a two-spool ducted fan engine and the other of the last 
stage which would be part of the fan turbine. Both designs assumed free- 
vortex flow, radially constant work, and 50 per cent stage reaction at mid- 
radius. Design I had a stage loading factor of 2.02, a flow factor of 0.75, 
and a pressure ratio of 1.9. A moderately high hub-to-tip diameter ratio 
was selected with which the rotor root reaction was limited to 15 per cent. 
The rotor exit angle was approximately constant at -62 degrees, the inlet 
angle varied from 10 to 55 degrees from casing to hub, and the rotor row 
relative velocity ratio varied from 0.4 to 0.75* The corresponding varia- 
tion in rotor total -pressure-loss coefficient varied from 0.07 at the tip 
to 0.39 at the hub. The mass flow weighted stage efficiency was computed 
to be 0.915 which compares well with the stage data of Figure 4. It is 
about one half of one per cent less than that predicted on a mean-line 
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basi s (see Fi g 15) . 

For Design II, a loading factor of 1.4, a stage flow factor of 
0.99, and a pressure ratio of 1.67 were selected. A lower hub-to-tip di- 
ameter ratio was also selected, and the resultant hub reaction was near 
impulse. The variation in rotor loss coefficient for the lower deflection 
blading varied from approximately 0.07 to 0.22. The mass flow weighted 
stage efficiency was found to be 0.898. This value is approximately 0.005 
less than the value obtained from either Figure 4 or Figure 15* 

Although it is unwise to base any firm conclusions on the results 
of a limited investigation, it would appear that the correlation can be ap- 
plied to elements of a stage design despite the fact that it was derived 
from a mean section performance analysis. 

Loss Factors 

While the correlation is to be made an integral part of a design 
analysis computer program, it is desirable to provide the program user with 
the facility to adjust the level of loss by input factors or to specify 
row loss coefficients as inputs. The correlation was based on achievable 
efficiency data, and the efficiency levels were adjusted to zero tip leak- 
age values. Thus, one of the principal corrections to be made for a tur- 
bine analysis will be for rotor tip clearance effects. The program will 
be supplied with the option to specify radial variations, of the loss fac- 
tor. These factors will be used to adjust the level of internally computed 

f 

loss coefficients. Therefore, I t wi 1 1 be possible t# 1/tc pease loss levels 

V. ^ ? 

in the vicinity of rotor tips to represent the additional TossNiue to tip 
clearance. However, the amount of the correction to be made and the radial 
extent of the additional loss is difficult to assess at this time. Various 
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empirical methods for the correction of stage efficiency for tip clearance 
effects exist, but none of these are specifically applicable to a stream- 
filament analysis. One of the possible problems associated with the stream- 
filament approach is that, if severe gradients of loss are specified, there 
will be no solution of the input specification. More specifically, in 
order to satisfy both radial equilibrium and continuity equations a re- 
verse flow, with negative values of meridional velocity, may be required. 

The manner in which radial distribution of additional losses and compensat- 
ing distributions of work output are specified, in order to maintain a 
mechanically acceptable standard of blade row geometry, can only be fully 
investigated when the computer program is available. 

Other factors which affect the level of loss are Reynolds number, 
blade row aspect ratio, and trailing edge thickness. As with tip clear- 
ance loss, various empirical correction methods exist for these factors. 
However, there are no generally accepted corrections; the reason for this 
is probably due to the fact that the detailed design of the blade is more 
important than the over-all interrow aerodynamics or the over-all geometry 
of the section against which these effects are often correlated. Similarly, 
the datum loss levels may be a function of the over-all design. For ex- 
ample, a constant section stator blade may have a distribution of loss which 
differs both in radial distribution and level from a conventional free-vortex 
stator. These differences may not be predictable on a simple stream-filament 
1^3 sis, since the complete distribution of the blade surface pressures within 
the blade passage may not be simply related to the blade row inlet and exit 
aerodynamics which are used as the basis for the loss correlation. 

While the design-point analysis may be undertaken to investigate 
over-all performance using the internal correlation of tota 1 -pressure-loss 
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coefficient and the input loss factors, it is believed that the program 
will have to be used in conjunction with test data from stages designed 
using the stream-filament method before suitable correlations for the ad- 
ditional loss factor can be developed. 


Kinetic-Energy-Loss Coefficients 

Although tota 1 -p ressure- 1 oss coefficients are used to express 
blade element performance in the analysis, there are, of course, other ways 
in which the row performance may be expressed. One of these alternatives 
is the kinetic-energy-loss coefficient^ , which is defined as follows: 


« i-(^Y 


( 60 ) 


where the subscripts denotes the isentroplc value of blade row exit ve- 
locity. The isentropic velocity is that velocity which would be attained 
in an isentropic expansion from the row inlet total pressure to the value 
of row exit static pressure. In Appendix I, it is shown that the total- 
pressure-loss coeff i ci ent , Y , and the kinetic-energy-loss coefficient,^ , 
are related by the expression 


Y ~ LLi 
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( 61 ) 


0 - 1 ) 

where is the tota 1 -to-stat i c pressure ratio at blade row exit. 

For rotor rows the relative total pressure is, of course, used in the 
pressure ratio. 

Since the solution of the flow conditions at each of the design 


planes will Involve the evaluation of the tota 1 -to-stat i c pressure ratio 
across the annulus, i t wi 11 be possible to readily accommodate kinetic- 
energy-loss coefficient as an optional design input. 
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DEVELOPMENT OF THE ANALYSIS PROCEDURE 


I n t roduct i on 

The solution of the flow field for the turbine design problem 
consists essentially of obtaining values of total pressure, total tempera- 
ture, and the two components of the absolute velocity at each selected 
streamline location in each of the design planes. Thus, for the analysis 
Po,To , VL , and are the principal variables. In the initial section 
of this report, all the equations which have to be satisfied in the analy- 
sis have been presented in terms of these variables. For an axisymmetric 
flow in an arbitrary annulus, the meridional slope and curvature of the 
stream-filament surfaces would have to be considered as analysis vari- 
ables, but in the current analysis these quantities are regarded as known 
and will be derived from the flow boundaries as defined by the specified 
annulus contours or directly specified as inputs. 

The solution of the flow field is obtained for the absolute 
values of P, , T c , ^ , and V M at each of the design stations. Conditions 
relative to the rotors at inlet to and exit from each rotating blade row 
are obtained using conventional turbine design techniques after the de- 
sign specifications and radial equilibrium have been satisfied in the ab- 
solute planes. Even though three types of design stations (namely, first 
stage stator inlet, stator exits, and stage exits) have to be considered, 
the design specifications for these stations differ, and various specifi- 
cation options are to be permitted for design-point analyses, the method 
of solution of the radial equilibrium and continuity equations has been 
developed so that the numerical procedures are identical for each design 
plane. 
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Before discussing the solution procedures, the specification of 
the inputs for a design-point analysis is reviewed. Following the presen- 
tation of the solution procedures, the outputs from the design-point an- 
alysis are also presented. 

Specification of the Design Requirements and Analysis Variables 

From the point of view of the numerical solution of the design 
problem, there is no real distinction between design requirements and 
analysis variables. Nevertheless, the following discussion subdivides 
the input into design requirements and analysis variables in the conven- 
tional manner. 


Design Requirements 

The type of unit, that is single or multispool, will be speci- 
fied by a simple indicator of the number of spools to be considered. The 
analysis variables will be read in on a spool -by-spool basis, and hence 
the i nd i ca tor % wi 1 1 principally serve to distinguish between new spool 
data and new analysis variables for the original spool. With this ap- 
proach the over-all program is considerably simplified, and the storage 
is relatively simply organized to accept up to eight stages on any spool 
and up to three spools in succession. It would have been possible to 
arrange to store data for more than one spool at any given time and thus 
permit successive analyses of alternative variables on multispool con- 
figurations. However, practical cons i derat i ons indicate that it would be 
uneconomical to attempt a detailed analysis of alternative design vari- 
ables for two or more spools. Preliminary analyses of a multispool con- 
figuration can be performed with any variations of the analysis variables 



considered as complete new cases. When an over-all design outline has 
been obtained, detailed analysis could then be performed on each spool in 
turn. The over-all procedure need not necessarily be slowed down by the 
need to establish the outlet flow conditions of a preceding spool to de- 
fine the inlet to any individual spool. The result of the preliminary 
multispool configuration can, of course, be used for more detailed in- 
vestigation of the analysis variables of second or third spools. 

The mass flow at machine inlet will be specified in Ibm per 
sec and in the case of a cooled turbine will be supplemented by the speci- 
fied cooling flow schedule. Coolant flows will be expressed as a fraction 
of the inlet mass flow and will be assumed to enter the main flow at the 
blade rows. Thus, if the coolant fractions for the first stator and 
first rotor are^^ y and.Xi^, respectively, then the continuity equation 
will be satisfied for>ur:, A*r T (t -f* , and (I **r CJ ) at design 

stations 0, 1, and 2, respectively. The temperature of the coolant at 
each point of admission will also be specified in general; an option will 
be provided so that these temperatures need not be specified if they are 
to be considered equal to those of the main stream. 

Since the input will be accepted spool by spool , the rotative 
speed and required power output of each spool will be simply specified 
by the rpm and horsepower output of each shaft. 

Inlet flow conditions of total pressure, total temperature, and 
flow angle will be specified against radius; the dimensions of these quan- 
tities will be 1 bf per sq in, deg R, degrees, and inches, respectively, 
in the computer program, parabolic interpolation of these quantities will 
be used. If only one value is specified, it will, of course, be assumed 
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constant; more complex distributions will require more data points to 
adequately describe them; if the inlet is that of a second or third 
spool, the distribution will be based on the output data from the pre- 
ceding spool. 


Analysis Variables 

The number of stages will be simply specified. This number 
will be used to identify the subsequent data for it will be assumed that 
each stage consists of a stator followed by a rotor row. 

The wall geometry will be specified by hub and casing radii at 
each design station. If these data are also to be used to compute wall 
slopes and curvatures, the axial spacing of these stations together with 
dummy stations ahead of the inlet and downstream of the final stage will 
also be specified. As an input option, streamline slopes and curvature 
in the meridional plane may be directly specified at selected radii. 

When slopes and curvatures are computed from specified annulus geometry, 
these quantities will be assumed to vary linearly between computed values 
at the annulus walls; when directly specified, streamline values of slope 
and curvature will be obtained by a linear interpolation or extrapolation 
of the specified data. 

The power output distribution will be specified both by stage 
and across the annulus for the individual stages. The stage-by-stage dis- 
tribution of output will be simply specified as a series of fractions of 
the total spool power. The distribution of power from hub to casing of 
individual stages will be specified using a nondimens ional power function 
which varies from zero at the hub to unity at the casing. These power 
functions will be specified for each streamline except for the special 
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case of uniform distribution of power output. The reasons for the selec- 
tion of this particular type of power distribution specification are dis- 
cussed fully later. 

Since the method of solution follows the streamline flow from 
the spool inlet, the basic specification of tangential velocities across 
the annulus will be made at stator rather than stage exits. The tangen- 
tial velocities at selected radial positions will be specified at each 
stator exit plane. However, since mechanically acceptable blade geome- 
tries are more readily obtained when the stator exit absolute flow angles 
are specified, the specification of tangential velocities or flow angles 

will be one of the program options. 

The row loss characteristics need not be specified if the in- 
ternal correlation of tota 1 -pressure-1 oss coefficient is used. However, 
to provide the flexibility for particular design analyses, the datum level 
of loss coefficients may be adjusted using loss factors in the input 
specification. These loss factors can be radius and row dependent. The 
use of the internal tota 1 -pressure- loss correlation and loss factors is 
undoubtedly the most useful option for a design analysis, the basic pur- 
poses of which are to determine the design geometry and predicted per- 
formance. However, additional options are provided. These are to specify 
(as functions of radius) tota 1 -pressure-loss coefficients, kinetic-energy- 
loss coefficients, and rotor or stage isentropic efficiency. The specifi- 
cation of isentropic efficiencies to represent the rotor row performance 
characteristics, in conjunction with a stator total-pressure-loss coeffi- 
cient, is an option which will not be frequently used for the analysis of 
new designs. However, the provision of this option will make the computer 
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program of greater value in that it could be used for the detailed inves- 
tigation of experimental data from which further information on the stream- 
filament loss characteristics of turbines will have to be obtained. 

Since the analysis will be performed with design station values 
of specific heat, a value of specific heat will have to be specified for 
each of the design stations. 

The analysis variables reviewed above complete the input speci- 
fication for the stream-filament analysis. However, when the effects of 
interfilament mixing are to be simulated in the analysis, the radius and 
row dependent mixing parameters will also have to be specified. 


Basic Equations and Fundamental Solution Technique 


In order to establish a procedure which is independent of the 
type of design station and selected input options, the analysis has been 
developed on the basis that the total pressure, the total temperature, the 
tangential velocity, and the meridional velocity are the unknowns, even 
though in certain circumstances one or more of these quantities will be 
directly specified. The basic equations which have to be satisfied are 
those concerned with 


1. Power output 

2. Radial equilibrium (radial momentum) 

3. The element performance (total pressure loss) 

4. Angular momentum - Euler work 


and 5. Mass fl ow conti nu i ty 

All the necessary equations, with the exception of the power output equa- 
tion, have been presented in the initial section of the report in terms of 
the analysis variables P c , 7^ , l/ R , and l/ fc . For the general design-point 
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analysis the variables P 0 , V^, and l# t are interdependent, but the fourth 
variable, 7^, is directly dependent on the specifications of the design 
variables. Thus, the evaluation of streamline values of total tempera- 
tures is regarded in the analysis as a preliminary to the solution of 
the problem proper and will be discussed prior to the method of solution 
for the three remaining variables. 

Evaluation of Total Temperatures 
The total temperature at the first design station, that is the 
first stator inlet plane, will be directly specified as a function of 
radius. A linear interpolation of these Input data will be used to es- 
tablish values of total temperature at the streamline positions used in 
the analysis. At the following station, a stator exit, the streamline 
values of total temperature will be unchanged from the corresponding 
streamline value at the preceding station for the purely streamline an- 
alysis, but will be modified in the manner previously discussed when in- 
terfi lament mixing or the addition of coolant have been specified. 

The total temperature distribution at the next station, a 
stage exit, will have to satisfy both the specified total power output 
and its distribution across the annulus. Since initially the distribu- 
tion of mass flow throughout the annulus is unknown until the distribu- 
tion of meridional velocity has been established, the power distribution 
will be specified by nondimens i ona 1 power functions versus the nondimen- 
sional mass flow function,/^) , which has been defined by Equation 21. 

If the total power output specified i s ffy (horsepower) , the 
total temperature drop AT 0 through the rotor must satisfy the equation, 

fP T - i&f A-r„ < 62 
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Normalizing Equation 62 with respect to the total power and the total 
mass floWj/M'j., leads to a definition of a power function ~1°(<*r(A)) which 
is expressed as 




5~sro 


I 


A T c c^.Wa.) 


where ^(a)\s the nondimensional mass flow function which varies from zero 
to unity between the hub and casing streamlines. Similarly, the power 
function will vary between zero and unity. Differentiating Equation 63 
with respect to the nondimensional mass flow function yields a general ex- 


pression for the total temperature drop; for the j t * 1 streamline, 

l Al ~) _ (eL-PfrrM)} 

^ 0y J AX ~~ r Crc f > l d /J 


(64) 


d Pm4A 

where aHa) ]S the local slope of the power function with re- 

spect to the nondimensional mass flow function. Thus, the power function 
versus mass flow function will be a basic specification for power distri- 
bution from which the total temperature drops are obtained. Hence,. the 
actual temperature distribution at stage exit will be obtained on a stream- 
line basis using the calculated temperature drop and the corresponding 
streamline total temperature at the preceding station. (Where mixing or 
the addition of coolant to the rotor row has been specified, the upstream 
values of total temperature will be modified as previously discussed.) 

The principal advantage of the use of this type of specification 
for power distribution, which in effect relates total temperature drops 
to streamline positions rather than radial positions, is that the power 
output equation (Equation 62) is automatically satisfied. It is there- 
fore possible to directly assign values of total temperatures to the 
streamlines without undertaking the iteration which would be necessary 
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if the specification involved power output as a direct function of radius. 
The type of specification is particularly well suited to a stream-filament 
analysis where the program user is concerned with the power output distri- 
bution between filaments rather than total temperature drop as a function 
of radius. This type of design will be characterized by the selected 
power function versus mass flow function relationship. If only 40 per 
cent of the total stage power output is required from the flow between 
the hub and mean streamline, one of the input data points will be 0.5, 

0.4; the slope of the power function (and hence the total drop along the 
hub contour) will be lower at the hub (0.0, 0.0) than at the casing (1.0, 
1.0). For a ’’constant work” distribution, the power function versus mass 
flow function will be linear between the points (0.0, 0.0) and (1.0, 1.0). 


Evaluation of Total Pressures and Velocity Components 
Except for the cases where the .total pressure and tangential 
velocities are directly specified or can be simply obtained from the 
Euler work equation, the total pressure and the tangential and meridio- 
nal components* of velocity are interdependent. The distribution of 
V and 14 at each design station must be such that they satisfy the re- 

r Jl ’ J)V 

quirements of radial equilibrium and continuity. The relevant equations, 
expressed in terms of the analysis variables, are repeated below. 
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The tangential velocity, as a function of radius, will be di- 
rectly or indirectly specified. At the first design station, which will 
be the first stage stator inlet of the design being analyzed, the tangen- 
tial velocity will be indirectly specified by the flow angle. Hence, at 
any radius 


Va* = V* 0 05 A* 


( 67 ) 


At the stator exit plane, when the option of specifying flow angle is se- 
lected, the relationship between V^j and v* I is, of course, similar to that 
of Equation 67. That is, 

Vwi - V m , A , (68) 


When stator exit tangential velocities are directly specified, then tan- 
gential velocity can be considered as a known function of radius. 

Vm = V W) M (69) 


Similarly, since streamline values of the total temperature drop across 
the rotor blade will be obtained from the specifications, the tangential 
velocity at stage exit can be also considered as a known function of radius 
if the radial locations of the streamlines are assumed. That is, 



(i. at; 




( 70 ) 


Considering the total pressure distributions, at the first de- 
sign station the total pressure will be specified versus radius. Hence, 
the relevant equation is 

Poo = Poo <*> <70 

At stator exit the total pressure must be obtained from the loss coefficient 
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equation, 


x, = 


Ppa Pc, 

Po, - P, 


( 72 ) 


Thus, the explicit expression for P a \ in terms of the analysis variables 
becomes 

R, = 


a 


Co 


i + y^D - 0 - 


(73) 


When the performance of rotor elements is specified by a total - 
pressu re~ 1 oss coefficient, the explicit equation for stage exit absolute 
total pressure, is derived from the loss coefficient definition which 

is repeated below 

Poxs ~ P pA- 


y* 




Hence, 


P - P 

' 02 ^ > O 
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(75) 

for / V)>' 

(Equation 

16), 


Using the previously presented relations 1 
Equation 75 can be reexpressed in terms of the analysis variables as 
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(76) 


where the isentropic value of rotor relative exit total pressure, P 0 xs , 
is given by Equation 15- The equations for stage exit total pressure are 
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considerably less complex when either of the alternative options of speci- 
fying stage isentropic or rotor isentropic efficiency is selected. These 


equations are obtained directly from Equations 12 and 13. That is, 


and 


Pox = Poo [ I - 

Po, * Po, [> - 


Toi ~?~o3. "I y-i 

Ijg I oc J 


(77) 
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T °° J (78) 


Since total temperatures can be obtained directly from the de- 
sign specifications, for these options the total pressures can be regarded 
as known functions of radius, Pz ’ Pox <*>. for the solution of the flow 
field. 

From the equations above it will be seen that in order to evalu- 
ate total pressures or tangential velocities it is, in general, necessary 
to first determine the meridional velocity. Since the distribution of 
meridional velocity must satisfy both radial equilibrium and continuity 
equations and these also involve the total pressure and tangential ve- 
locity, the mc^t obvious method of solution would be the solution of three 
simultaneous equations. However, the radial equilibrium equation is a 
differential equation involving the derivatives UV Yd<., and 

k dP °W . Hence, the two additional equations must be obtained from 
the differentiation with respect to radius of the tangential velocity and 
total pressure equations. Since the three equations can only be solved 
for the derivatives, the actual values of the variables must be obtained 
from the simultaneous solution of the continuity equation. 

Before discussing the numerical techniques which are used for 
the solution, the following section presents the particular forms of the 
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differential equations which are used in the solution. 


The Differential Equations 

In presenting the equations, it will be assumed that the design 
stations to be considered are the first stator inlet, the first stator 
exit, and the first stage exit. The following stations in a multistage 
machine can, of course, be considered as merely repeating the problems 
presented by the first stator exit and first stage exit planes. 

The radial equilibrium equation (Equation 65 ) is of the form 


dVm 


I JLPo 


d V*. _ 


( 79 ) 


where C„ , C lt , C, 3 , and C ( * are coefficients which can be assigned values 
at each point in a particular design plane once a value of meridional ve- 
locity has been selected. The variables 7^, tosf\, -A*and A, and the con- 
stants & , , and T are assumed to be known quantities. 

The differentiation of the appropriate equation for total pres- 
sure (which is either Equation 71 , 73 , 76 , 77 , or 78 , depending on the 
station being considered and the selected option) leads to a differential 
equation also of the form 
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( 80 ) 


where the coefficients C u , , C l3 , and C ^ are again in terms of quan- 

tities which are readily obtained from the design specifications for an 
assumed value of the meridional velocity. 

Similarly, differentiation of either of the alternative expres- 
sions for V L , (obtained from Equations 67 , 68 , 69 , or 70 ) will produce 
differential equations of the form 


6k 


o ^ ^ S. dJl ^ Vu 

5/ ~C&b- ^35. P, + ^33 'cDtT ~ & 2*t (8i) 

The coefficient C lA i s j n each case zero. 

Although the actual differentiations involved to determine the 
coefficients are not in themselves complex, the expressions for the co- 
efficients are in some cases quite complex. For convenience of presenta- 
tion, the actual coefficients are given in Appendix II. 


Technique for Solution 

While the coefficients of the three differential equations will 
differ depending on the station being considered and the selected input 
option, once these coefficients have been evaluated a unique technique can 
be used for the solution. Essentially, the problem becomes one of obtain- 
ing a meridional velocity distribution which simultaneously satisfies the 
radial equilibrium and continuity equations. Selecting an initial value 
of meridional velocity at one streamline position, the local values of the 

coefficients of the set of equations can be obtained. These equations are 

(JiV* 

then solved for the derivative . Then, using standard forward differ- 

ence techniques, the value of meridional velocity at an adjacent stream- 
line is obtained using the originally selected meridional velocity and the 

, . dvi 

derivative • Using the new value of meridional velocity, consistent 

values of total pressure, tangential velocity, and the coefficients are 

obtained for this streamline using the appropriate equations previously 

cl V ^ 

presented. Thus, the derivative can be obtained at this new stream- 

line also. The over-all process is repeated until the meridional velocity, 
total pressure, and tangential velocity have been determined at each of the 
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streamlines used in the analysis. Using Equation 18 , the continuity equa 
tion, the mass flow for the distribution is evaluated. Since the distri- 
butions will have been based on an assumed value of meridional velocity at 
one point in the flow field, the mass flow computed will in general differ 
from the specified value. Hence, the assumed value of meridional velocity 
will have to be modified iteratively until the starting value is consis- 
tent with the continuity requirement. 

Over-All Solution Procedure 

The over-all design analysis proceeds from known inlet condi- 
tions station by station through the turbine. The basic calculations are 
performed using grid points within the flow field which are defined by an 
even number of equal-flow stream filaments. Since initially the flow dis- 
tribution is unknown, the initial streamline positions are estimated from 
equal areas for each filament. Hence, streamline positions have to be re- 
located after each solution of radial equilibrium and continuity until a 
converged solution for streamline positions has been obtained. Included 
in this major iterative loop will be an iteration on streamline values of 
total-pressure-loss coefficients when the optional specification of kinetic- 
energy-loss coefficients has been selected. 

When the meridional velocity satisfies the radial equilibrium 
equation, the specified design variables, and the continuity equation 
within a preset tolerance, new streamline positions, and where applicable, 
loss coefficients are obtained. Revised values of the streamline depen- 
dent variables required for the coefficients of the three differential equa- 
tion are then obtained. The solution of the radial equilibrium and conti- 
nuity equations is then repeated until the streamline positions, and loss 
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coefficients if ki net i c-energy- loss coefficients are specified, have con- 
verged to within preset tolerances. 

Having obtained the basic solution at one design station, 
streamline values of all the relevant aerodynamic parameters in both 
relative and absolute reference systems are readily obtainable using 
conventional turbine design procedures. Among the quantities computed 
will be those necessary as input for the solution of the flow field at 
the following design station. These will include, where applicable, re- 
vised streamline values of total pressure and total temperature when the 
addition of coolant flow and/or interfilament mixing has been specified 
for the downstream blade row. 

From the point of view of a numerical solution, the following 
design stations are solved in an identical manner. The only basic dif- 
ferences between stations and input options are in the evaluation of the 
streamline coefficients of the set of three differential equations and 
in the selection of the initial estimate of the meridional velocity. 

Because it is possible to have two solutions to compressible 
flow problems, it is advisable to commence the simultaneous solution of 
the radial equilibrium and continuity equations at a streamline which is 
most representative of the flow in the annulus. Hence, a mean stream- 
line is selected, which equally divides the flow in the annulus. In prac- 
tice this selection complicates the logic of the computer program in that 
the solution of the meridional velocity distribution has to proceed to 
each of the two boundary streamlines in turn. Nevertheless, for stator 
exit planes in particular, the variation in absolute Mach number across 
the annulus will be sufficiently large that convergence of the required 
solution will be best achieved when the meridional velocity is reestimated 
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at the most representative streamline for the flow field. 

When the flow angle is specified, for example at the first sta- 
tor inlet or stator exit planes, both subsonic and supersonic solutions 
are possible. At stator inlet it will be assumed that only the subsonic 
solution is of interest, and the initial meridional velocity will be se- 
lected to correspond to a Mach number of 0.4. At stator exit planes, it 
will be necessary to specify which of the two solutions is required. If 

the subsonic is chosen, the initial estimate of meridional velocity will 
be based on a mean Mach number of 0.8; for supersonic solutions the start- 
ing point of the flow iteration will be a Mach number of 1.2. 

When tangential velocities are specified at stator exit and 
when they are indirectly specified as is the case at stage exit planes, 
two solutions are again possible. However, only one is of real interest, 
since the second will correspond to a design in which the axial component 
of velocity is supersonic. For these cases the first estimate of meridio 
nal velocity will be based on a stator exit angle of 60 degrees or a 
rotor relative exit angle of -60 degrees. For all design analysis of 
practical interest, the numerical solution will converge to that for 
which the axial component of Mach number is subsonic even though the 
absolute Mach number may be either subsonic or supersonic. 


Resul ts of the Analysis 


Since on completion of the basic solution, streamline values of 
total pressure, total temperature, the velocity components, and the flow 


angles will have been obtained at all design stations, the computer pro- 
gram output can be arranged to have any of the turbine design parameters 
as output. In line with good practice, the output will, of course, 
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contain a print out of all the input specifications of design requirements 
and analysis variables. Velocity triangle data will be given at each 
streamline station by station. These data include: 

Radi us of s treaml i ne 

Meridional, axial, and tangential velocity components 
Absolute velocity and flow angle 
Blade speed and relative flow angle. 

The state of the gas at each of the streamline positions will be defined 
by: 

Absolute total pressure and total temperature 
Rotor relative total pressure and total temperature 
Absolute and relative Mach numbers 
Static pressures and static temperature. 

Where interfilament mixing and/or the addition of coolant has been speci- 
fied, values of absolute and relative total pressures and total tempera- 
tures which are assumed to exist at the inlet to the following row for 
the purpose of calculating the following design station flow field will 
also be tabulated. 

Following each complete stage a performance summary for the in- 
dividual streamlines will be presented. These data will comprise; 

Stator and rotor velocity ratios as indicators of the section 
react i ons 

Stator and rotor tota 1 -pressure- 1 oss coefficients 

Stator and rotor blade row efficiencies (defined as where 

€ is the kinetic-energy-loss coefficient 

Rotor and stage isentropic efficiencies. 

The stage output data will be completed by a tabulation of mean values 

of: 
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Stator and rotor blade row efficiencies 
Stage work output in Btu per Ibm 

Stage tota 1 -to-tota 1 and total -to-static isentropic efficiency 
Stage blade-to-jet speed velocity ratio. 

Except for the last item which will use the mean streamline values, the 
above items will be computed using mass flow weighted values of pressures 
and temperatures. 

At the conclusion of the design analysis for a spool, mass flow 
averaged quantities will be presented for the multistage unit. These 

output data will include: 

Spool work and power 

Over-all total -to-total and tota 1 -to-stat i c pressure ratios 

Over-all total -to-total and tota 1 -to-stat i c isentropic effi- 
c i enc i es. 

An over-all blade-to-jet speed velocity ratio based on the over-all spool 
pressure ratio and a mean blade speed will also be presented. 

For*the purely stream-filament analysis, that is in the ab- 
sence of specified mixing or a coolant flow schedule, all the output 
quantities will be obtained using standard turbine analysis formulas. 

When mixing or coolant flows are specified, the definitions of stream- 
line work and efficiency and mass flow weighted mean values of efficiency 
will be as defined in the earlier section of the report concerned with 
the modification of the basic analysis procedures. Since complete data 
on all the total temperatures and total pressures will be available as 
output, alternatively defined stream-filament efficiencies which would be 
consistent with the total -pressure-loss coefficients used in the analysis 
in the case of specified mixing within blade rows will be readily 
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calculated from the output. In addition, mass flow weighted efficiencies 
for cooled turbines can be redefined and reevaluated if an alternative 
definition of efficiency is preferred. 


71 


REFERENCES 


1. Smith, S. F. , "A Simple Correlation of Turbine Efficiency", J. Royal 
Aero, Soc . , vol. 69 , July, 1965* 

2. Stabe, R. G. , et al , Cold-Air Performance Evaluation of a Scale-Model 
Fuel Pump Turbine fo r the M-l Hydrogen-Oxygen Rocket Engine (NASA TN 
D-3819), National Aeronautics and Space Admi ni st rat i on , Washington, 

0. C. , February, 1967. 

3. Stewart, Warner L. , et al, "A Study of Boundary-Layer Characteristics 
of Turbomachine Blade Rows and Their Relation to Over-All Blade Loss , 
Trans. ASME, The American Society of Mechanical Engineers, Series D, 
September, I960, p. 588. 


72 


NOMENCLATURE 


Symbols 

Description 

Units 

A 

Angle of streamline slope in 
the meridional plane 

deg 


Specific heat at constant 
pressure 

Btu/lbm deg R 

-P 

Total blade surface diffu- 



s i on 

-- 

£ 

Kinetic-energy-loss coeffi- 
cient (= |- ) 



»/ 


£ 

Alternative kineti c -energy - 
loss coefficient for compres- 



s i b 1 e flow ( = 1 - (ft Vp/ ( ^ v ^) 

-- 

tf 

Arbi t ra ry f unct i on 

— 


Constant in Newton's law 

Ibm/lbf ft/sec^ 

k 

Entha 1 py 

Btu/1 bm 


Total power in horsepower 

hp 

J~ 

Mechanical equivalent of heat 

ft lbf/Btu 


Index on streamlines 

— 

/TV 

Number of streamlines 

-- 

P 

Total pressure 

psi 

p 

Static pressure 

psi 

r 

Nondi mens i ona 1 power function 



_ MT r Tc h . 

- * T ‘ 

— 

P 

Stage power output 

Btu/ sec 

R 

Stage reaction 

-- 

R 

Gas constant 

ft Ibf/lbm deg R 


Radi us 

i n or f t 

5 

Pi tch 

i n or f t 
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Description 


Units 


Symbols 

X 

** 

it 

v 

w 

yCtr 

/Ur(x) 


y 

fi 

y 

$ 

*l R 

f 

<f> 

f 

Sl 

Subscri pts 

c 


Z. 1 — 

Total temperature 

Trailing edge thickness 


Blade speed 
Veloc i ty 

Work extraction along a 
st reamtube 

Hass flow rate 

Cooling flow fraction of to- 
tal. annulus flow 


Nond i mens i ona 1 mass flow func- 
tion Jx '- 
Ur 


’ A*r r ft 


_ ( f±\\ - { V M U>sA A 

Jr. L 1 SgeTcpTc J 


deg R 
in or ft 
ft/sec 
ft/sec 

Btu/lbm 

Ibm/sec 


Mixing parameter 

Tot a 1 -pressure- 1 oss coef f i c i ent 

Flow angle 

Ratio of specific heats 
Row deflection 
Rotor isentropic efficiency 
Stage isentropic efficiency 
Densi ty 

Stage flow function 

(- v*/& ) 

Stage loading function 

(= /e ) 

Rotational speed 

Descri ption 

Coolant 


deg 


deg 


lbm/ft 3 


rpm 
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Subscr ? pts 


Descript i on 


CL. 

Casing 


Exi t 


Hub 

/CVL. 

Inlet 

YK 

Meri di onal 

N 

Stator ("nozzle 11 ) 

R 

Rotor 

A' 

Radi al 

$ 

Stage 

S 

1 sent ropi c 

T 

Total 

U. 

Tangential 

X 

Axial 

o 

Stage inlet 

1 

Stator exit/rotor inlet 


Stage exi t 

Superscri pts 

Descr i pt i on 

/ 

Relative to rotor 

* 

Subsequent to mixing 


Subsequent to coolant addi 


Mean or mass flow weighted 
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APPENDIX I 


THE RELATIONSHIP BETWEEN TOTAL-PRESSURE-LOSS 
COEFFICIENT AND KINETIC-ENERGY-LOSS COEFFICIENT 


The losses In a blade row, or an element of a blade row, can be 
expressed in a number of ways. In the current analysis NREC has adopted 
a total-pressure-loss coefficient rather than a kinetic-energy-loss co- 
efficient, which is preferred by some turbine designers. There are no 
significant difficulties associated with the introduction of both types 
of loss coefficient into the turbine design analysis program, since the 
two coefficients are simply related for a given value of blade row exit 
Mach number. In the following, the relationships are derived and a 
third coefficient, which is a true kinetic-energy coefficient for com- 
pressible fluids, is also briefly discussed. 

The tota 1 -pressure- 1 oss coefficient is defined as follows: 

PcZS ~ Po 


y = 


on 


( 1 - 1 ) 


Pan - 

where K*s >s .the isentropic value of total pressure at blade row exit, 
P OJL i 5 the actual total pressure, and Pn is the row exit static pressure. 
For stator rows the total pressures are absolute total pressures, and the 
isentropic total pressure at exit is equal to the row inlet total pres- 
sure ^ ( i . e. , fez*, — Pot )• For rotor blades, the total pressures are 

relative values, and the isentropic total pressure at row exit will only 
equal the inlet total when there is no change of radius between inlet and 

exi t. In general , 


Pc 


oz 5 


+ 


Po, 


i sl 




( 1 - 2 ) 
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(1-3) 


The definition of the kinetic-energy-loss coefficient is as follows: 

<L ~ I- 


v; 




where is the velocity at blade row exit and V;^ is the isentropic ve- 
locity at blade row exit, which is defined as the velocity which would be 
obtained for an isentropic expansion from the inlet total pressure to the 
actual static pressure at row exit. As with the total-pressure-loss co- 
efficient, absolute quantities are used for a stator row and relative 
values for a rotor. 

It should be noted that the coefficient defined by Equation 1-3 
is not strictly a kinetic-energy coefficient for compressible fluids, 
since the actual density at row exit will not equal that obtained by an 
isentropic expansion to the same row exit static pressure. This point is 
discussed later. 

The relationship between £ and ^can be expressed in many forms, 
but since the total-pressure-loss coefficient involves total and static 
pressures, a logical choice for the relationship is in terms of total-to- 
static pressure ratio at row exit. 

Equation 1-1 can be rewritten as follows: 



Hence, the ratio of i sentropic- to-actual total pressure at row exit is 
given by 


/s2S 

PoA. 


i +■ y (i- 



(1-5) 


Equation 1-3 can be expressed in terms of tota 1 -to-stat i c pressure ratio 
as 
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( 1 - 6 ) 



Hence , 


JH 

7 


_ _ (MhiM 

e " ' -\tr 


(1-7) 


Equation 1-7 can more conveniently be expressed as 


e = 






y-\ 

v 


( 1 - 8 ) 


i 


Substitution of the expression for total pressure ratio given in 
tion 1-5 into Equation 1-8 produces the following expression for 
terms of Y and . 


Equa- 
^ i n 



[i + y( 

1 -£)] 

r-i 

1 

7 

1 +- 7( 


RSe) 

7-1 


<o ^ 

Equation 1-9 can, of course, be rearranged to obtain an explicit 
sion for / in terms of £ and Pc2 /a- is, 


(1-9) 

expres- 


y = 


fe) -(MZtf'T*- < 

i- F - 


( 1 - 10 ) 


Po 


The expressions given in Equations 1-9 and 1-10 can be written in terms 
of exit Mach number, critical velocity ratio, or total-to-static tempera- 
ture using the standard thermodynamic formulas. Regardless of the par- 
ticular form of the relationship, i t wi 1 1 be seen that total -pressure-loss 
coefficients may be readily transformed to a kinetic-energy coefficient 
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and vice versa for any given value of tota 1 -to-stat ic pressure ratio (or 
Mach number). With tas an input to a design program, or contained in 
the program as a basic correlation, values of the total -pressure-loss 
coefficient used in the solution of the radial equilibrium equation (which 
forms the heart of the analysis procedure) can be obtained in an iterative 
procedure. This iterative procedure will essentially involve obtaining 
row exit total pressures from previously assumed or computed row exit 
flow conditions. It should be pointed out that this iteration is simi- 
lar to that which is required when a correlation for total pressure loss 
is used in the analysis program. 

It can be shown from Equations 1-9 and 1-10 that if Y is as- 
sumed independent of Mach number, then ^ wi 1 1 decrease with Increasing 
row exit Mach number, and conversely if is assumed independent of exit 
Mach number, V will increase with Mach number. In practice, which of 
the two coefficients is more nearly constant with varying Mach number is 
not established. In the context of a turbine design-point analysis, it 
is implicitly assumed that any turbine design requirement will be suit- 
ably bladed, and hence the behavior of loss coefficients with Mach number 
determined from tests of given blade geometries is not strictly relevant. 
For example, a blade designed for low Mach number operation may exhibit 
increasing loss coefficients with increasing Mach number, whereas a sec- 
tion specifically designed for high Mach number operation may have the 
reverse loss coefficient versus Mach number character i st i c. 

In the earlier discussion of the kinetic-energy coefficient, 
it was stated that the expression given as Equation 1-3 assumes that the 
isentropic density,^, is equal to the actual density, The actual 
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kinetic-energy-loss coefficient should therefore be written as follows for 


compressible flow; 

// / _ f* 

£ ^ 1 TTvg o-io 


Hence, assuming the exit static pressure and exit 
equal for both isentropic and actual expansions, 
pressed as follows: 


t 


// 


vC /ft* 

Vzl V P -3 


total temperature are 
Equation 1-11 can be ex- 


1H 

I r (1-12) 


Therefore, from Equations 1-3 and 1-12, 


u 


(1-13) 


ii 


C - I — U ~ \ Pozs 

From Equations 1-7 and |-13 the relationship between & and C. in its si 
plest form becomes 


m- 


c 7/ = e 


Pcrz \ y 


= 


(1-14) 


If it is assumed that €. is independent of exit Mach number, will de- 
crease with Mach number, whileY will increase with Mach number. This 
third coefficient is of interest therefore in that it does provide an 
intermediate between the differing schools of thought, the one consider- 
ing tota 1 -pressure- 1 oss coefficient independent of Mach number and the 
other which assumes C is independent of Mach number. 
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APPENDIX I I 


COEFFICIENT S FOR THE EVALUATION OF 
MERIDIONAL VELOCITY DISTRIBUTION AT ANY DESIGN PLANE 

The meridional velocity distribution must satisfy the radial 
equilibrium equation (Equation 65). Since this equation involves deriva- 
tives of , and Vu_ , the method of solution selected is the simul- 

taneous solution of three differential equations in which tK* V J_ 

11/ cCr Po 

3,1 cCr are un k nowns ’ ^he two additional equations are obtained from 

the differentiation of the appropriate equations for total pressure and 
tangential velocity. If a value of is assumed at an initial stream- 
line, the solution of the set of three equations will yield a value of the 

l , / TL 

derivative . Hence, using a forward-difference technique, the value 

a ^ an ac U'acent point in the flow field may be evaluated and the 
solution repeated. Having obtained a complete solution for the annulus, 
the mass flow passing through the design station is computed. If the con- 
tinuity equation (Equation 66) is not satisfied, the initial value of me- 
ridional velocity is reestimated and annulus flow conditions recomputed. 
Thus, the radial equilibrium and continuity equations are satisfied using 
an iterative procedure. 

To simplify the logic of the computer program, a standard proce- 
dure is adopted for the solution of the flow field at each design station. 
The different types of design stations and the various optional specifica- 
tions are taken into account by modifications to the twelve coefficients 
which appear in the three differential equations (Equations 79, 80 , and 81 ) . 
This appendix presents these coefficients. The three types of design sta- 
tions are considered in turn; where there are optional specifications, 
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these alternative forms of the coefficients are given. 

When the second differential equation is derived from the defi- 
nition of the total-pressure-loss coefficient (at stator and stage exits), 
the expressions for the coefficients C 1( , Cx\ , and Ci* contain addi- 
tional coefficients C Y , , Cy» , and C Y<r . The values of these additional 
coefficients will depend on the selected correlation of total-pressure-loss 
coefficient. The actual expressions for Cy, , Cy-j , and for the se- 

lected loss correlation are presented in Appendix III. 


First Stator Inlet 

The coefficients of Equation 79 are 


Cir- l*o 








_ _ avC + Lv*cV\CY_l 4^ 

-y- T "let 

1 m 


The total pressure will be a specified function of radius. Hence, the co 
efficients of Equation 80 are 


C-oj - 0 
Cm.* i o 

c a - o 



tion of flow angle A c with radius. Hence, from the differentiation of 
Equation 67, the coefficients of Equation 81 are 
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'3i 


/~a v\ cca /)c - 


3.V 


1*16 


C 32 . = CO 

I ^ 


^ l (A — taviA {Lt'n fkclflA 

W/k c ; U i cTFj 

Stator Exi t 


The coefficients of Equation 79 are 

= l 0 

Ca jV a ,\ V^^TC^j 

C/3, ^ aV-i 

C, ( Vui%- ^, X )jz- - 

r Y7v\ ” »Oi CU * 

The coefficients of Equation 80 are obtained from the differentiation of 
Equation 73. In the analysis it is assumed that the derivative d Y* can 
always be expressed as 

dY, = Cv. clUO +Os.lVL + Cr . 

to 3; a? * 


Thus, the coefficients of Equation 80 are 

0 ,= (M(1l Y„ + (SVi- £lc r , 

■a^RTc, U<A rj 


Cv, - | 0 


c, . OMIT V. + 1& 

,% RTli 



Po 




’9* 


r C c + ^“ ,lr ^ ^ - 1 !ioV i -i?\ c Y 

Im \ Pc A Pc) + 
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If the tangential velocity is a specified function of radius, the coeffi- 
cients of Equation 81 are 

C 3t . o 

O4 x. ■* o 
- VO 

Q* - 

(,LT 

Alternatively, if the flow angle, f\ , is the specified quantity, the co- 
efficients of Equation 81 are obtained from the differentiation of Equa- 
tion 68. Thus, the coefficients are 
Cjj - 

X Vmi 


O 

x l o 



Stage Exi t 

The coefficients of Equation 79 are 


G „ -- l-o 

ClX ^ * ^ t T( t Tc "} 


c>$ = aV* jl 

...Itoflx - + (vC»v.„0i- ^ 

it- cvr 


c It , . K 


V, 


p\ 


r 


If either the option to specify rotor isentropic efficiencies or to specify 
stator isentropic efficiencies as a function of radius are used, the local 
total pressures can be calculated from Equation 77 or 78. Thus, the total 
pressure can be regarded as a known quantity and the coefficients of Equa- 
tion 80 are 


Cxi = 0 




I O 
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C ^ 

. Poi 7CF 

If the performance of the rotor is specified by a loss coefficient, the 
coefficients of Equation 80 are derived from the differentiation of an 
expression for fcjL • Equation JG can be reexpressed as 

Pco. — Pci A I 

U — 

% * 


(i i-D 


where 


and 


3i 


I i + * lZ 1 



x 

«■-! 


^ ! ' 1 ~ = / V-'_ | ( 4- U^-<3u x V (czl ^Tj 

f/ v ' L a.i -TClollx 


S' ■ t ■ H-T- 


TCjp IC>J 

t — \ /ux l llM 

L- *1 iy j J 




Hence from Equation I 1-1 

J_ . Po^ _ __L . 

Pci- cif df cW 


•+- 


t d > - 1± ^ w 


+ \ X - - 1 / 1 - 9 }\ d Yr 
^ M ,V/ X 1 


Thus, when the differentials are evaluated 
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_L. cl fc— _ _L.il [ft/ 4 . J f - (uJ--a, x ) j_ 

P oi . ctf foi X<* 2‘y^,| oi . L T C i c'It J 


~ i ii Ml 

2 ^ ft To/ 

- [ av,v 


2 . u z ft — 2 a/.clVji- - 2\ r ^z 

M 


ft_ — (u./~ jl<t A ui)_l_ . t 
)cz 


3* R ‘ 


-{sLi+Vut)^- ‘IIL- 

ftz cXr 


J clC , 

L 37 ^ 

- I fki )( l ~ CVi ■ ^ +- ^Tj a Vuz. -f- Cy ^-1 

UxA fez a. ax cw J 

where it is assumed that the derivative tlYa. can be expressed as follows; 

tX-T 

‘■1 ^ s. Cxi i t ^ f- Gyj 1 1 Vu j. 4. CYia. 

nr ar tc? 

Thus the coefficients of Equation 80 will be 


Cm, y 


+ 


' 2^4 f.Tz 


IfciX 1 - *)■ 


C»= 1-0 


g, , v4MH - -4j-€ 


^ft 1 


Soft 1 ' 




£ _ J_ . Jft' , J ( c -?» x il -3 u, m T- - (u^-uAL 4 <o/| 

^ Pc/ ctr ‘i<j t RX/ l Tci ttf J 


- i 

ftW t 


a - 2 V Ul n- - (u4- ^Vu^-L Jj« 

lei cl^f 


f ^ ( gJ (^MC)X. cljoz -/&JV i-it-Yv* 

£3. RT X T «. ^ ^ Pcz^ XX/ 
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Since the streamline total temperature can be readi ly calculated 
from power output function, the radial variation of tangential velocity, 
can be calculated from Equation 70. Thus, S/u-l is in effect a known func- 
tion of radius and hence the coefficients of Equation 8l are 


Cv * G 

t-Vl. ~ o 


C& = lo 



dr' 
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APPENDIX I I I 


THE COEFFICIENTS OF THE 
TOTAL-PRESSURE-LOSS COEFFICIENT DERIVATIVE 

The second differential equation used in the simultaneous solu- 
tion of iV* is obtained from an expression for the local value of total 

aC^- 

pressure. When the design station is a stator or stage exit and the loss 
of total pressure through the preceding blade row is computed from a 
tota 1 -pressure-1 oss coefficient, this second equation will involve the 
derivative dlY-v or tlY«l . In Appendix II it has been assumed that these 

■acr p 

can be expressed as follows: 

— (- v f| d 4- 4- Cy<*- ( I 1 I “ 1 ) 

7CF P 

where the coefficients Cyj , Cyj. and W**- will appear in the expressions 
for C* , Cxs » and Cxq. . 

If the tota 1 -pressure-loss coefficient is a specified function 
of radius, or the over-all solution is an iterative one in the case of the 
specified kinetic-energy-loss coefficient option, these coefficients are 
as fol lows : 

Cy, = Cyj a O and Cy 4 — d YCv ) 

tit" 

When the local value of the total -pressure-loss coefficient is at its maxi 
mum allowable value of Qtj , then 

Cy, _ Cy^ =■ Cyy.= 0 when Y - a<^ 

When the tota 1 -pressure- 1 oss coefficient is obtained from an internal cor- 
relation, the expressions for Cy, , Cy^ » an( * will depend on the 
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particular correlation being used. For the computer program the loss co- 
efficient will be expressed as follows: 


Hence, 


Y- i, k <3 

ft 


dX = Y j 4" dLfc +. -t ^4^ 

cW~ 


(I I 1-2) 


f Y' ^ -h-Ldh-S-jU 

( J| cU- fa ZC7 ■$> ~dCr j 


(III -3) 


The values of Qy f , (y^ , and are then obtained by equating the right- 

hand sides of Equations MI-1 and IM-3. 


Stator Exi t 


The loss coefficient is defined by 


^ X/ - 4 (~ ^2Jb. - Y Z a ^1 when -> ( 1 I I -4a) 

( 0-1+ + Os' Ccs/3, ) 


Y n _ j (- + qi (*)** 1 

+- ayO 


when ^^<*3 (Ml-4b) 


where ^ , is the additional loss factor which can be specified as a function 
of radius. The form of the individual derivatives will depend on whether 
the flow angle or the whirl velocity is the specified quantity 


and the value of 

V, 


Flow Angle Specified and 




C Y = - X<<K x Vc 
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C\ 5- ;\ { co/i 


£\'a. 



t ^ 'iw% r f\ % cc? /l, ti^L‘ 

;jrv cCr 


Stage Exi t 

The definition of rotor loss coefficients is essentially identi- 
cal with that used for the stators (Equation I I I -4a and I I I -4b) ; relative 
quantities and the relevant design station indices must, of course, be 
substituted throughout* Note that ^>2 is by definition, a negative 
number 


U/J ~ ii ( h'.vy-s, - hw ) [a , t- a z [ _ a ^j| / 

1 7- L — =-r if ¥*>43 (III -5a) 

(^^q. as J V/ 

Y ft , (m . B) 

V'/ ' 




Jill 


Q; 


c 


rj - Y* 


r ~ a^' _ 
L V'- 



-x, 




+ 
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Stage Loading i-actor 



Stage Flow Factor, Vtc/lL 
Note: Reproduced from Reference 1 
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FIGURE 9 - REDUCED LOSS COEFFICIENTS WITH 
ADDITIONAL EXIT ANGLE CORRECTION VERSUS VELOCITY RATIO 
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Stage Isentropic Efficiency 










Stage I sentropi c 



Stage Flow Factor, 'V/./lL- 


FIGURE 11 - A COMPARISON OF TEST DATA EFFICIENCIES WITH PREDICTION 
VALUES USING ALTERNATIVE LOSS COEFFICIENT CORRELATIONS = 1 . 5 ) 
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Experimental Data (From Fig 4 ) 

Predicted Using Correlation A 
(Equation 58 a) 

Predicted Using Correlation B 
(Equation 58b) 

Predicted Using Correlation C 
(Equation 58c) 
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Stage Isentropic Efficiency 
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FIGURE 13 - A COMPARISON OF TEST DATA EFFICIENCI ES WITH PREDICTION 
VALUES USING ALTERNATIVE LOSS COEFFICIENT CORRELATIONS fa.TCpA.T./to =2 . 5) 
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Stage Loading Factor, 0© ^ 



FIGURE 14 - PREDICTED EFFICIENCY CONTOURS 
.I^raWrTSS COEFFICIENT 
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Stage Loading Factor, j3* 
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Stage Flow Factor, V x / 


FIGURE 15 ^ PREDICTED EFFICIENCY CONTOURS 
BASED ON ROW LOSS COEFFICIENT CORRELATION B 
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